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Abstract. In previous work, we showed that the solution of certain systems of discrete inte- 
grable equations, notably Q and T-systems, is given in terms of partition functions of positively 
weighted paths, thereby proving the positive Laurent phenomenon of Fomin and Zelevinsky for 
these cases. This method of solution is amenable to generalization to non-commutative weighted 
paths. Under certain circumstances, these describe solutions of discrete evolution equations in 
non-commutative variables: Examples are the corresponding quantum cluster algebras [3], the 
Kontsevich evolution '10' and the T-systems themselves [S]. In this paper, we formulate cer- 
tain non-commutative integrable evolutions by considering paths with non-commutative weights, 
together with an evolution of the weights that reduces to cluster algebra mutations in the com- 
mutative limit. The general weights are expressed as Laurent monomials of quasi-determinants of 
path partition functions, allowing for a non-commutative version of the positive Laurent phenom- 
enon. We apply this construction to the known systems, and obtain Laurent positivity results for 
their solutions in terms of initial data. 
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1. Introduction 

1.1. Discrete integrable systems and positivity. Discrete integrable systems are evolution 
equations in a discrete time variable, which possess sufficiently many conserved quantities (discrete 
integrals of motion) so that they can be solved explicitly in terms of some discrete Cauchy initial 
data. 

Cluster algebra mutations p!5] can be viewed as discrete evolution equations for certain com- 
mutative variables, where at each step, the evolution is a rational transformation. The evolution 
is known to result in variables which are Laurent polynomials as functions of the initial data, even 
though the evolution equations are rational transformations |16j . 

Cluster algebras are algebras over commutative variables (the cluster variables) in which sets 
of variables (the clusters) are associated with the nodes of a regular n-tree, and cluster variables 
on connected nodes are related by an evolution called a mutation. One can choose a cluster at 
any fixed node to be the initial data for the evolution. Fomin and Zelevinsky conjectured that in 
terms of any initial data, the cluster variables at any node are positive Laurent polynomials [15], 
This is known as the positivity conjecture. 

The Ar Q-systems and T-systems are discrete integrable evolution equations which are a subset 
of a cluster algebra structure [21^16]. These systems were solved, and thus their positivity property 
proved, in ^ (for the T-systems, see [5j for generic boundary conditions). The approach of [8] 
and [9] involves the interpretation of the solutions of the systems as partition functions for paths 
on some discrete target space, with weights explicitly related to the (discrete Cauchy) initial data 
of the evolution equations. 

The T-system can be considered as a particular non- commutative version of the Q-system, in 
which we adjoin a non-commuting element to the commutative algebra generated by Q-system 
solutions |9] . The T-system equations are matrix elements of the cluster algebra mutations twisted 
by the non-commutative element. The solutions are then expressed in terms of partition functions 
for paths as for the Q-system, but with non-commutative weights. 

In all these examples, the positive Laurent property follows from the fact that the weights are 
explicit positive Laurent monomials of the initial data. 

The aim of the present paper is to formalize this underlying path structure in the more general 
case of arbitrary non-commutative weights, with the idea of defining non-commutative (integrable) 
cluster algebra mutations with a built-in non-commutative positive Laurent property. 

We note that there exists an example of a non-commutative cluster algebra with these proper- 
ties: the quantum cluster algebra [3], where the cluster variables g-commute. This is the mildest 
form of non-commutativity, and in this case, much more can be said about the structure of the 
solutions. 
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1.2. The Aj. T-system and Q-system. The T-system is a system of discrete evolution 
equations for the commuting variables {Tij^k}- It can be written in the fohowing formQ 

(1-1) 7i,j,fc+iTij,fc-i = Ti^j+i,kTi^j-i,k + Ti+ij^kTi-ij^k, 'i^ < i < r,j,k e Z, 

with the convention that Tqjj^ = Tj._|_i = 1. 

Equation p.ip (with a sign change on the right-hand side) was first introduced [2^, in 
the context of the generalized Heisenberg spin chain in statistical mechanics. With appropriate 
boundary conditions, it is the "fusion relation" satisfied by its transfer matrices. Its solutions are 
the g-characters of Uq{5ir+i) pTj [27]. 

A closely related system is the Ar Q-system [22], 

(1-2) Ri^k+iRi,k-i = Ri,kRi,k + Ri+i,kRi-i,k, 1 < i < r,k e Z, 

with i?o,fe = Rr+i,k = 1- This system is obtained from the T-system (II. Ij) by dropping the index 
j, for example, by considering solutions of (II. ip which are 2-periodic in j, Tjj+2,fc = 21,j,A; for all 
i,j,k. Then Ri^k = ?i,i+fcmod2,fc is a solution to ([O])- 

Remark 1.1. Given the appropriate boundary conditions (Rifl = (— 1)LV2J j ffi^g system is satisfied 
by the characters Qi^n = {—'^)^'^^'^^Ri,n (n > 0) of the irreducible Ar-module with special highest 
weights which are multiples of one of the fundamental weights. Using Equation ()1.2p . these are 
expressed as polynomials in the characters of the fundamental modules, identified as {Qi^i : 1 < 
i < r}. 



In reference to the origin of the equation, we refer to Equation (11. 2p as the A^ Q-system 
throughout this paper, without imposing this special boundary condition. 

To solve such a system means to express the variables for any n in terms of the initial data set. 
For example, valid initial data sets for the Q-system are sets of 2r variables {Qi^mi, Qi,mi+i ■ 1 ^ 
i < r} such that |mj — mj+i[ < 1. For fixed values of {i,k) or {i,j,k), Equations (ll.lj) . (ll.2j) are 
mutations of cluster variables in a cluster algebra. We will therefore investigate the systems from 
the point of view of local mutations of their initial data. 

1.3. Example: Commutative and non-commutative Ai Q-systems. The simplest case of 
the Q-system occurs for r = 1. This example was previously considered from a different point of 
view by several authors [H [281 ISSl [l3]. Let us recall the solutions of both the commutative and 
non-commutative versions of the Ai Q-system [71 [10] . We wish to stress the similarity of their 
path solutions, which relies on the existence of a sufficient number of conserved quantities and a 
linear recursion relation involving these. 

Consider the recursion relation for the commutative variables {Rn}, 

(1.3) Rn+lRn-l = (Rnf + 1 

The system has a manifest translational invariance. We thus compute the solution n > in 
terms of the initial data {Ro,Ri). 

The system is integrable: The quantity c = ^j;^^ ~^ r r R^+i independent of n. Moreover, 
the solutions of (|1.3p obey the linear recursion relation Rn+i — cRn + Rn-i = 0. Let 

yi = Ri/Ro, y2 = l/{RoRi), Vs = l/vi = Ro/Ri- 
and define the generating function 

Mn 1 



(1.4) Fo{t) = Y,i 



n=0 ^0 1-17%^ 



-'^The relative minus sign which appears in the original version [2] [24] is normalized away in this paper. 
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This is the partition function for weighted paths on the integer segment [1,4], from 1 to 1, and 
with steps ±1, with weights 1 for step i i + 1 and tyi for step i + 1 — )• i = 1, 2, 3. As mentioned 
above, this interpretation makes positivity manifest, as the y^'s are exphcit Laurent monomials 
of the initial data {Rq , Ri). 

The non-commutative Ai Q-system is in a class of wall crossing formulas introduced by Kont- 
sevich, and involves the variables R„ € a skew field of rational fractions in the variables 
R^^jR]*^^. It can be written as |10j 

(1.5) R„+i(R„) ^R„_i = R„ + (R„) ^ 

The goal is to express R^, n > in terms of the initial data (Ro,Ri). 

The system is integrable, in the sense that there are two conserved quantities [10], independent 
of n: 

C = (Rn+l) "'"RnRn+l(Rn) ^, K = R„+i(R„) + (R„+i) "^(R„) + (R„+i) "^R„ 

and the solutions to (II. 3p obey the linear recursion relation Rn+i — KR,„ + CR„„i = 0. Define 
yi = Ri(Ro)~\ y2 = (Ri)"HRo)"^ and = (Ri)~^Ro, with ysyi = C. Then the generating 
function 

CO _-|^ 

(1.6) Fo(t) = rR„(Ro)-^ = (l - (1 - (1 - ty^r'ty2)~%i) • 

n=0 

This can be interpreted as the partition function of weighted paths on the integer segment [1,4], 
from 1 to 1, and with steps ±1, with non-commutative weights 1 for i — )• i + 1 and tyi for i + 1 — )• i, 
i = 1,2, 3. In this case, the weight of a path is the product of step weights, ordered according to 
the order in which they are taken. 

We see that the non-commutative Ai Q-system, despite its complexity, has essentially the same 
solution as the commutative one: it is expressed in terms of paths on the same graph, with similar 
but non-commutative weights related to the initial data, and with the condition that the order in 
which weights are multiplied is that in which the corresponding steps are taken. 

We recall that in the commutative case, the mutation of initial data in terms of which i?„ is 
expressed may be implemented via local rearrangements of the finite continued fraction expression 
(II. 4p for Fq (t) [7] . This property was shown to carry over in the non-commutative case |10j , where 
mutations of the initial data are obtained via local non-commutative rearrangements of the non- 
commutative finite continued fraction expression (II. 6p for Fo(t). 

1.4. Plan of the paper. The example above suggests the following strategy for exploration of 
the non-commutative evolutions. Instead of starting from an evolution equation, we start from 
the "solution", namely some path models on the same target graphs as those considered in the 
solution of commutative systems, but with non-commutative weights. It turns out that the par- 
tition functions of these paths have a structure similar to that observed in the commutative case. 
In particular, we will define non-commutative mutations connecting them via non-commutative 
finite continued fraction rearrangements. Moreover, by investigating the relation between the 
coefficients of these continued fractions (i.e. the weights of the path models) and the coeffi- 
cients of their series expansion, we will be led naturally to express results in terms of the quasi- 
determinants introduced in |18) . Our main result is the construction of a chain of path models 
with non-commutative weights, all related via mutations, and whose partition functions involve 
a basic set of non-commuting variables, in a manifestly positive Laurent polynomial way. 
The paper is organized as follows. 

We first review in Section [2] the solution of the (commutative) Ar Q-system, in a new path for- 
mulation related to that of [7] (the equivalence is described in detail in Appendix [X]). Each initial 
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data, and therefore each graph and associated set of weights is coded by a Motzkin path of length 
r — 1. The corresponding partition functions of weighted paths display a simple finite continued 
fraction structure, and the mutations are shown to be expressed through local rearrangements. 

Section[3]is devoted to the non-commutative generalization of this structure, and involves paths 
on the same graphs but with non-commutative weights. Mutations are shown to be expressed 
again as some non-commutative local rearrangements of the finite continued fraction expressions 
for the partition functions of paths, giving rise to simple evolution equations for the corresponding 
non- commutative weights. We then define non-commutative variables by discrete quasi- Wronskian 
determinants of the coefficients of the series expansion of the partition functions. These variables 
are shown to satisfy an evolution equation (the discrete non-commutative Hirota equation. The- 
orem I3.27P playing the role of the Q- or T-system in the non-commutative setting. Finally, we 
express the weights of all the mutated path models in terms of these discrete quasi- Wronskian 
variables (Theorem 13. 33p . as manifestly positive Laurent monomials. 

Sections|3]and[5]are applications respectively to the T-system, as an example of non-commutative 
Q-system, and to the quantum Q-system, in relation with the corresponding quantum cluster al- 
gebra [3]. In both cases, we compute the discrete quasi- Wronskian variables in terms of the initial 
data, and obtain Laurent positivity theorems for the solutions of the corresponding systems. 

In Section 6, we give some further examples of non-commutative evolutions which have a Lau- 
rent positivity property. Section 7 considers the problem of generalizing of the Lindstrom-Gessel- 
Viennot theorem to the non-commutative setting. We conclude with some remarks about the 
possible generalization of cluster algebra relations and variables to the generic non-commutative 
setting. 

Acknowledgments. We would like to thank S. Fomin for discussions, and L. Faddeev for 
pointing out Refs.[T2| I13j to us. We thank the Mathematisches Forschungslnstitut Oberwolfach, 
Germany, Research in Pairs program (Aug. 2009) during which this work was initiated. PDF also 
thanks S. Fomin for hospitality at the Dept. of Mathematics of the University of Michigan, Ann 
Arbor (Spring 2010). RK thanks the Institut Henri Poincare, Paris, France, for hospitality during 
the semester "Statistical Physics, Combinatorics and Probability" (Fall 2009), and the Institut de 
Physique Theorique du CEA Saclay, France. PDF received partial support from the ANR Grant 
GranMa, the ENIGMA research training network MRTN-CT-2004-5652, and the ESF program 
MISGAM. RK is supported by NSF grant DMS-0802511. 

2. Continued fraction solutions to Q-systems 

2.1. The Ar- Q-system. We consider the set of discrete evolution equations 
(2.1) Ri^k+iRi,k-i = Ri,k + Ri-i,kRi+i,k, ^ < i < r, k e Z 
with the boundary conditions i?o,fc = Rr+i,k = 1 for all A; G Z. 

2.2. Initial data sets. For fixed a and n, the variable Ra,n+i is a rational function of i?/3,n with 
|a — /3| < 1 and Ra,n-i- Therefore the recursion relation (|2.ip has a solution for a given set of 
initial data, which must include the 2r variables -R(m) = {Ra,mc ^ Ra,mci+i '■ 1 < a < r} where the 
sequence m = (mi, m^) is such that \mi — mi+i| < 1. That is, initial data sets of the Q-system 
are in bijection with Motzkin paths of length r — 1. 

Any variable Ra,k which is a solution of (|2.ip can then be expressed as a function of the initial 
data set R{ni). Moreover, due to translational invariance of Eq. (j2.ip . the expression for Ra,n as 
a function of R{m) is the same as that of R^^n+k as a function of R{k + m), with the translated 
Motzkin path A; + m = (A; + mi, ...,k + m^). We may therefore restrict ourselves to a fundamental 
domain of Motzkin paths modulo translations, M = {m|Min(mj) = 0}. 
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In [21], we showed that the relations ()2.ip are mutations in a coefficient-free cluster algebra, 
defined from the seed cluster variables {xq,B), where 



where C is the Cartan matrix of the Lie algebra Ar- This seed corresponds to the fundamental 
Motzkin path, m = (0, ...,0). 

Not all mutations in the cluster algebra are of the form ()2.ip . However, all the variables 
{Ra,n '■ I < a < r, n £ 7j} are contained in the subalgebra made up of the union of seeds 
parametrized by the set of Motzkin paths of length r — 1. Thus, we consider the corresponding 
sub-cluster algebra, restricting to the nodes of the Cayley tree that are labeled by Motzkin paths, 
and to the mutations that connect them. 

We note that mutations in the sub-cluster algebra under consideration here have the following 
property: fii is given by one of the Q-system equations if and only if it maps ii(m) to i?(m'), 
where m'j = rrij it (5ij. That is, it maps a Motzkin path to another Motzkin path which differs 
from it in only one node. 

2.3. Path partition functions. 

2.3.1. Positivity from weighted paths on graphs. In [7j, we first eliminated from the Q-system 
(12. ip the variables Ri^k i > 1, A; € Z, in terms of the Ri^m n £ 7^, hy noting that the discrete 
Wronskians 



actually satisfy the same recursion relation (12. ip . as a consequence of the Desnanot-Jacobi identity 
for the minors of any given matrix (also a particular case of the Pliicker relations, see (7j for 
details). As moreover VFo,n = 1 and Wi^n = Ri,n, we conclude that Wi^n = Ri,n for all i > 0, 
n G Z. 

The Q-system therefore boils down to a single non-linear recursion relation for Ri^m in the 
form Wr+i^n = 1- In Ref.[7j, we explicitly solved this equation, by using its integrability. The 
latter implied the existence of a linear recursion relation for the Ri^^s, with constant coefficients 
(independent of n, and expressible explicitly in terms of initial data). Exploiting this, we further 
showed that Ri^n could be related to the partition functions of certain weighted path models on 
target graphs Fm (see definition in Appendix lA]l . indexed by Motzkin paths m of length r — 1, 
the weights being directly expressed as positive Laurent monomials of the corresponding initial 
data i?(m). 

Finally, the quantity Ri^n, expressed as the determinant Wi^n i\2.2\i was shown to also be related 
to path partition functions, namely that of collections of i strongly non- intersecting paths on the 
same target graphs Fm, with the same weights. This followed from a slight generalization of the 
Lindstrom-Gessel-Viennot theorem |25] |20j . 

Moreover, we obtained an explicit expression for the generating function of the variables \^R\^n '• 
n > 0} as the partition function of paths on the graph Fm, in the form of an explicit finite 
continued fraction. 

In this section, we present a similar but simpler description in terms of paths on another, more 
compact graph Gm, with related weights. The equivalence of the two formulations is proved in 
Appendix Rl 




(2.2) 



Wi^n — det (i?l,n+i+l-a-6 



l<a,fe<i 
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2.3.2. The weighted graphs G^. For fixed r, consider the weighted graph G^. The graph itself 
depends only on r. It has r + 1 vertices labeled 1, ...,r + 1 and oriented edges connecting vertex 
i to itself, to i — 1 and to i + 1. This is pictured below, 



Gr 



Q 



r + 1 



1 



where each edge connecting two different vertices is to be considered as a doubly oriented edge. 
The map m H- Fm i-> Gm is described in Appendix lAl 

To make the contact between the Gm path formulation and that on FjYj, we need a few defini- 
tions. We start from a collection of 2r + 1 formal variables, y(m) = (yi(m), y2r+i (m))- These 
will eventually be identified with the skeleton weights of the graph Fm (see Appendix [A|) . 

Definition 2.1. To y(m) we associate bijectively the collection y{m) = (yi(m), y2r+i(in)). 
This is defined according to the local structure of the Motzkin path m as follows: 



(2.3) 



m = ty2i{ty2i+ir^+'-"'', 

ty2i^i-:j^, mi+i=mi-l; 
y2i-i = \ ty2i-i+ty2i, mi+i=mi + l; 

ty2i-i, otherwise. 



Example 2.2. The three basic examples are as follows: 

• //m = mo = (0, 0), the fundamental Motzkin path, we have jji = tyi, i = 1,2, ...,2r + l. 

• If m = mi = (0, 1,2, ...,r — 1), the strictly ascending Motzkin path, we have 

hi-i = t{y2i-i + y2i), hi = i^2/2i2/2i+i(l <i <r -1), 

with xja = tya for a >2r — 1. 

• If m = m2 = (r — 1, r — 2, 0), the strictly descending Motzkin path, 



m-i = ty2i-i 



with ya = tya for a > 2r — 1 . 



y2i 

y2i+i 



y2i 



y2i 

y2i+i ' 



(1 < i < r - 1), 



We now attach weights to the edges of the graph Gm- The edge connecting vertex i to vertex 
j carries a weight Wij := Wij{ni), where 



with yi{m) as in (|2.3 



2.3.3. Partition functions of paths on Gm o,nd continued fractions. Each path p on the graph Gm 
from node i to node j carries a weight wt{p) which is the product of the weights Wi^j of the edges 



8 PHILIPPE DI FRANCESCO AND RINAT KEDEM 

i j traversed by the path. The partition function of paths from i to j is the sum of the weights 
of ah paths p € "PijiGm) which start at vertex i and end at vertex j: 

Zij{G^) = J2 wt{p). 

Once expressed in terms of {i/i}, the weights of Definition 12.11 introduce terms with negative 
signs. The following lemma shows why this does not introduce minus signs in the generating 
function of paths. 

Lemma 2.3. The partition function Zij{Gm.) o.s a power series in t, has coefficients which are 
positive Laurent polynomials in {yi\i- 

Proof. Let I± = {i : rrii = rui-i ± 1}. The partition function for paths on Gm is the same as that 
on the graph G'^, where each loop at vertex i such that i € /_ is decomposed into two loops as 
follows: 



b 

a — b 



b 



The corresponding piece of is depicted on the right. In this case, o = ty2i-i and b = 
y2i/y2i+i- We call the two corresponding loops on the right-hand side if (with weight a) and 
(with weight —b). Note that the weight of i'^ is equal to — tfj+i^j, as shown in the picture. 

We define an involution on paths p G J'ij (G^) as follows. 

• v{p) = p p does not pass through any (i G /_) nor contains any sequence of steps 
i -> i + 1 i (i G /_); 

• Otherwise, consider s, which is either the first step of p passing through some , or the 
first sequence of the form i i + 1 i [i ^ /_) inp, whichever comes first. The path (p{p) 
is obtained from p by replacing s by the other sequence (a step through £~byz— + 7>z 
and vice versa). 

Clearly, if is an involution on J't,j(G"jn). Its effect is to reverse the sign of the weight, wt{ip{p)) = 
—wt{p) if (p{p) ^ p- There is therefore a pairwise cancellation of contributions of paths such that 
p 7^ 'p{p)- The terms which are left in the partition function are the fixed points p = <f{p) of the 
involution. These have no step through any loop (i G /_) nor any sequence i i + 1 i 

{iel-). 

The remaining weights are all positive Laurent polynomials of the y^'s, and the lemma follows. 

□ 

Given the graph Gm with weights Wjj(m), the partition function of paths from vertex 1 to 
itself can be written as a finite continued fraction as follows. 

Definition 2.4. The Jacobi-type (finite) continued fraction is defined inductively as 
(2.4) jM(xi,...,X2.+i) = - -, J«(^' ^ 



l-Xi-X2J(^ 1)(X3, ...,X2r+l) ' 1 
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Lemma 2.5. The generating funetion of weighted paths on from and to the vertex 1 can be 
written as 

(2.5) Zi,i(G^) = jM(y(m)). 

Proof. By induction. □ 

We regard J^'"^(y(m)) as a power series in t. As we have seen in Lemma 12.31 the coefficients of 
t"" are positive Laurent polynomials in y(m). 

Alternatively, the positivity of the generating function Zi^i(Gm) may be recovered by a simple 
rearrangement of the continued fraction j('')(y(m)). Starting with the bottom term in the con- 
tinued fraction and moving upwards, one iteratively rearranges the fraction for each i G /_ using 
the identity 

b 7. b + ^^u 

(2.6) a-- + - ^ = a + 



c I — c — u I — c — u 

with a = ty2i~i, b = ty2i and c = ty2i+i- (Here, u denotes the already rearranged remainder of the 
fraction). This eliminates any negative weights in the fraction, at the cost of a modification of the 
branching structure of the fraction (the fraction u appears in two places), and the introduction 
of new, positive weights coming from the term ^-u. 

The fraction J^'^\y{m)) can be re-expressed in a canonical, manifestly positive form as follows. 
Going though the fraction iteratively from bottom to top, whenever i I^, apply the identity 
(j2.6p . Otherwise, if = rui-i + 1, apply the identity 

/ r be b 
(2.7) a + b+- = a + 



l-c-u ' l-j^ 

with a = ty2i-i, b = ty2i and c = ty2i+i- (Here, u is the already rearranged remainder of the 
fraction.) 

The resulting rearranged form of J^^')(y(m)) has the property that all its coefficients are Laurent 
monomials of y(m) with coefficient t. 
We call this form canonical. 

Definition 2.6. The canonical form of the fraction J^^^(y(m)) is the fraction resulting from the 
use of the rearrangements (j2.6|) and (j2.7p . such that the final fraction has coefficients which are 
all Laurent monomials of y{m), with coefficient t. 

Another important generating function is the partition function of paths from l' to l' on the 
graph Lm of Appendix [Al It can be written as follows: 

(2.8) F™(i) = l + tyi(m)jM(y(m)). 

Remark 2.7. The canonical form of F^{t) itself (with the identity (12. TP applied also at the top 
level, with a = 0, 6 = 1 and c = yi(m)j is related directly to the path interpretation of Fm{t) 
as generating weighted paths on the graph (see Appendix [A], Eq AA.2\ ): while the yi's are the 
skeleton weights o/Fm, the redundant weights have been produced by the rearrangements ()2.6p . 
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Example 2.8. The fraction corresponding to the Motzkin path (0, 1,0) for r = 3 reads: 

yi 



^0,1,0 (t) 



l + t 



l-t{yi + y2)-t^- 



y2yi 



1-t 



yi 



i-t- 



■l-t!/5-*T^' 



Note the occurrence of a new "redundant" weight ty^yQ/y^. 

2.4. Mutations. Without loss of generality we now restrict ourselves to forward mutations of 
m, : m 1-^ m', with m'^ = rnj + 6j^i. We also restrict ourselves to two special cases: 

• Case (i): fii : m^, mj+i) = (a, a, a + 1) i-^ (a, a + 1, a + 1) 

• Case (ii): /Uj : (mj_i, mj, mj_i) = (a, a, a) (a, a + 1, a). 

These two cases are sufficient to obtain any Motzkin path, up to a global translation by an overall 
integer, from the fundamental Motzkin path mo = (0,0, ...,0). 
We also define the action of a mutation on Z(y): 



(2.9) 



(2.10) 



Case (i) 



Case (ii) 



m'. Define y' : 


= y{m') as a function of y := y{m) as follows: 




'y2i~i 


= 2/21-1 + y2i 


■ 1 


y'2i 


= y2i+iy2i/{y2i-i + y2i) 






= y2i+iy2i-i/ {y2i-i + y2i) 




'y'2i-i 


= 2/24-1 + y2i 


) : < 


y'2i 


= y2i+iy2i/{y2i-i + y2i) 




y2i+l 


= 2/2i+l2/2i-l/ (2/2i-l + 2/2i) 




,2/21+2 


= 2/2i+22/2i-l/ (2/2i-l + 2/2i) 



while all other weights are left unchanged, 

We then have the natural definitions of y(m') in terms of y(m'), using Equations (12. Ih . 
Theorem 2.10. Under the mutation /ij, the function j('")(y(m)) is transformed as 



(2.11) 



jM(y(m)) 



jM(y(m')), i>l; 
l + tyi(m')j(''Hy(m'), i = l. 



This has an immediate corollary. 
Corollary 2.11. Under the same mutation fii, F^nit) is transformed as: 

(2.12) F^it) = 



Fn,'{t), i>l; 
l + tyi{in)F^,{t), i = l. 

To prove the theorem, we use two rearrangement lemmas in addition to Equations (|2.6|) and 



Lemma 2.12. For all a,b we have 
(2.13) 1 + 



1-a-b 1-^ 
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Lemma 2.13. For all a, 6, c, d, with a + h invertible, we have 

h a' 

(2.14) a + 

where a',b',c' are defined by 

(2.15) a' = a + b 



1 ^ 1 _ b' 

i-d J- l-c'-d 



be 



a + b 

We now turn to the proof of Theorem 12.101 
Proof. Let Ji(m) be defined inductively as follows: 

y2i 



ac 
a + b 



Jj(m) = y2i-i + 



, (i>0), Jr.+i(m) = y2r+i(m), 



1 - Ji+i(m) 

with the convention that y_i = and = 1. Note that Jo(m) = j'^'^\y{^))- We refer below to 
the "level" i of the fraction to denote the part of the fraction J^**) where Jj appears. Let 

a. = ty2i-z{m), 6 = t7/2j-2(m), c = ty2i-i{m), d = ty2i{m), 
e = ty2»+i (m) , / = ty2i+2 (m) , u= ■ 

together with their primed counterparts, which are functions of m'. 

We now consider the effects of the mutations of Def. 12.91 on j('')(y(m)). 

• Case (i): We first transform the Jacobi-like form into the local canonical form for Jj_i(m) 
by applying (12. 7p at the level i, namely: 

Ji_i(m) = a + ^ — = a + 



l-c-d 



de 
1-e-fu 



1 



l-fu 



Next, we apply Lemma 12.131 to the second level of this fraction, namely: 

b 



Ji_i(m) = a + 



1 



with c' = c + d, d' = de/c' , and e' = ce/c' 



l-e'-fu 

Finally, we put this back into Jacobi-like form, by applying (|2.7p backwards at the level 
i — 1. The complete transformation reads: 

b . , . be' 



Jj_i(m) = a + 



l-c-d 



Ji_i(m') = a + b' 



l-e-fu 



1 



l-e'-fu 



i + 1 eQt 



f 



i + 1 

de 
i 

I 

i - 1 




Case (ii): We apply directly Lemma 12.131 to the level i, namely: 

b b 



Ji_i(m) = a + 



1 



a + 



l~e~fu 



1 



-, with < 



a'H j—u 

e 

l-e'-f'u 



(c' 


= c + d 


d' 


= de/d 
= ce/c' 




[f 


= cf/c' 
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We bring this back to Jacobi-like form by applying first (j2.6p backward on its first level, 
and then (I2.7P backward on its level zero, namely: 

b . be' 



Ji_i(m) = a + 



1 



a + b + 



1+ 



d' 



e' l-e'-f'u 

The complete transformation reads: 

b 



Ji_i(m) = a + 



1 -c 



Ji-i{m') = a + b + 



e' l-e'-f'u 



be' 



1-e-fu 



^ + e' l-e'-f'u 



i + l 



i - 1 



O 



i + l 



i - 1 



i + l 






/' 


et J 


H + l 




e' 




► i 

<.bc' 







Note that for z > 1, we always have Jo(m') = Jo(m), or equivalently -Fni(t) = F^i{t). When 
z = 1, we have a = 0, 6 = 1 and in all cases e' = yi(m'), so the final formula amounts to 
Jo(m) = 1 + t Jo(m')yi(m'), or equivalently -Fm(i) = 1 + tF^i {t)yi{in) . 

This proves the theorem and its corollary. □ 

2.4.1. Explicit expression for the weights yi{m). We can now identify the definitions of the mu- 
tations on the initial data variables fii : R{m) i-^ R{m') and the action on weights fii : y{m) — > 
2/(m'). 

Corollarv 12 . 1 1 I mav be used as a definition of the series Fm{t), by using induction over mutations, 
once we specify the initial condition. We can do so by giving the value of (t) as a function 
of i?(mo), the variables corresponding to the fundamental Motzkin path. This will define the 
initial values yj(mo) of the weights. Using Definition 12.91 we can then determine y(m) for any m 
inductively. 

In Ref . j7j , we used the integrability of the Q-system (|2.ip to show that the generating function 

-fo(i) = t"'Ri,n/Ri,o 

n>0 

is expressed explicitly in terms of the variables R{mo) as the finite continued fraction: 

yi 



Fo{t) = l + t 



l-tyi- t- 



y2 



ty-i-t- 



a2r+l 

where 

Ri,lRi-l,0 Ri+l.oRi-1,1 

(2.16) y2i-i = -—^ y2i = 

Comparing this expression with that for F^it) (j2.8p . one sees that the choice of weights yj(mo 
yi as in (|2.16p gives 

Fr^,it) = Foit) 

We can now derive explicit expressions for the variables y{m) for any m using Definition [27 
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Theorem 2.14. For any Motzkin path m = (mi, 1712, ...,m.r), the weights yi{m) defined by ()2.9p - 
()2.10p with the initial conditions yi(mo) as in (j2.16p are expressed in terms of the variables R{m) 
as: 

(2.17) y2i~l(m) = '^''"''"^^ , y2i(m) = di,rn,,+l-^-^,^+l(-^-^,^-l)~^ 



where 



and 




— — it rui = TUi+e + £; 
1, otherwise. 



T -^j+l,m^i— l,m— 1 



Proof. By induction under mutations, from mo to any Motzkin path m. We exphcitly check the 
relations p.9p - ()2.10p on the expressions ()2.17p . by use of the Q-system. □ 

Corollary 2.15. The series F^{t) is the following generating function for solutions of the Q- 
system: 

(2.18) Fm(t) = j;m /Ri ,mi 

n>0 



Proof. Whenever a mutation //j : m 1— t- m', i > 1 is apphed, the function -Fm(i) = Rm'{t) 
is preserved, hence its series coefficient remain the same. When i = 1 however, we have the 
transformation Fm(t) = 1 + tyi{m)F^i{t) = 1 + t ^"^^^^ F^i {t) , and therefore the coefficients of 
the series -Fm'(*) change according to (j2.18p . □ 



3. NON-COMMUTATIVE WEIGHTS AND EVOLUTIONS 

In this section, we introduce a non-commutative version of the path partition functions and 
weight evolutions which appear in the Q-system. Most of the key ideas of the previous section 
generalize to the non-commutative setting because paths are themselves non-commutative objects, 
and respect ordering of non-commutative weights. 

3.1. Paths with non-commutative weights and non-commutative continued fractions. 

Let {yf 1 : 1 < i < 2r + 1} be non-commutative variables and let A be the skew field of rational 
functions in {y^^, y^^,^^^}. Let Gm denote the graph used in Section r2.3.2[ We associate to 
each edge of a non-commutative weight, which is an element in A[t]. 

A path on a graph with non-commutative weights is defined to have a weight equal to the 
product of weights of the edges traversed by the path, written in the order in which the steps are 
taken, from left to right. For example, if the path is of the form p = (ii, 22, ■■■,in), then its weight 
is wt(p) = Wj^^j2Wj2,i3 • • • where Wjj- is the weight of edge from vertex i to vertex j. 

3.1.1. Non-commutative weights. For each Motzkin path in the fundamental domain, m € M, 
let y(m) = (yi(m), y2r+i(m)) S be a collection of 2r -|- 1 non-commutative vari- 

ables. Let I± = {i : rrii = mj_i it 1} as before. To y(m) we associate the variables y(m) = 
(yi(m), ...,y2r+i(in)) G defined as follows: 

y2i-i = t{y2i-i + y2i), 92i = i^y2i+iy2i, if i g /+; 
(3.1) y2i-i = iy2i-i - y^i+iy2i, y2j = y^,+iy2i, ifiG/_; 

y2i-i = iy2i-i, y2j = iy2i, otherwise. 

Here, we have omitted the reference to the Motzkin path m. Note that these are non-commutative 
versions of the variables y defined in Equation (j2.3p . 
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Finally, we define the weight on the edge from i to j of to be Wj j, where 

= y2i-i(m), Wi,i+i = 1, Wi+i,i = y2i(m). 

3.1.2. Non-commutative partition functions. Let 'J'a,b be the set of paths on Gm from vertex a to 
vertex b. Each path p E [Pq f, has a non-commutative weight 'wt{p) G ^1 defined like in Section 
12.3.31 as the product of step weights, in the order in which the steps are taken. The partition 
function, or generating function, of paths from a to 6 is defined to be 

(3.2) Z{G^)a,b= Yl ^*(P)- 

P&a,b 

3.1.3. Non-commutative continued fractions. 

Definition 3.1. For any x = (xi, X2r+i) with Xj G A[t], we define the non- commutative 
Jacohi-like finite continued fraction J^''^(x) inductively by J^"^^(xi) = (1 — xi)~"^ and 

(3.3) j(^)(xi,X2,...,X2fc+i) = (l-Xi - j('=~l)(x3,X4,...,X2fc+i)x2) \ 

As in the commutative case, we may consider such a fraction as an element of ^[[i]], that is, 
as a power series in t. 

It is clear that Z(m)i^i = j('")(y(m)) in terms of the weights (|3.ip . We also define the generating 
function Fm(t), in analogy with the commutative case, to be 

(3.4) F^(t) = l+tjM(yi(m),...,y2r+i(m))yi(m) 

Lemma 3.2. The coefficient oft"" in the power series expansion ofF^{t) is a Laurent polynomial 
in y(m) with non-negative integer coefficients. 

Proof. The proof is identical to the proof in the commutative case (c.f. Lemma [2.3p . since weights 
respect the order of steps in the paths. Suppose i G /_, so that Wj^j has a negative term. 
Decompose the loop i ^ i with weight Wjj = y2i-i into two separate loops from i to i, one, if , 
with weight ty2i-i and the other, /~, with weight — y^j^;^y2i = Using the involution ip 

defined in Lemma 12.31 on paths on the resulting graph G'^ , we again conclude that we have a 
cancellation of all paths p ^ ^{p) in the partition function. □ 

An alternative proof of Lemma l3.2l uses the structure of non-commutative continued fractions. 
One may simply rearrange the continued fraction expression (13. 4p (from bottom to top) whenever 
2 G /_, by use the following non-commutative version of the identity (12. 6p . For any a, b, c, u G A: 

(3.5) a - c^^b + (1 - c - u)"^c"^b = a + (1 - c - u)"^(b + uc^^b) 

Once again, this is done at the cost of introducing multiple branching of the new fraction, since 
u now appears in two places. Moreover, new "redundant" weights appear, from the new term 
uc^^b. These are, however, positive Laurent monomials in y(m). 

For use below, we note that we also have an analogue of the identity (12. 7p : For any a, b, c, u G A, 

(3.6) a + b + (l-c - u)"^cb = a + (l - (1 - u)"^c) "^b 
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3.1.4. Canonical form of continued fractions. As in the commutative case, the fraction Fm(t) can 
be rearranged into a "canonical" form, a (finite, multiply branching) continued fraction whose 
coefficients are all Laurent monomials in y, with coefficients equal to t. 

This is done as follows. From the bottom to the top level of the fraction (z = r, r — 1, 1), 
we rearrange the fraction whenever rui ^ mi-\. Let (a, b,c) = t(y2i-i, y2i) y2i+i) and use (|3.5p 
if i G /_, (|3.6p if i € /+, to rearrange the fraction. This includes a last "level zero" step of the 
continued fraction as well (case i = 0), which corresponds to the conventions y_i = and yo = 1. 

The resulting "canonical" continued fraction is interpreted as the generating function for paths 
from and to the root vertex on the graph Tm (see Appendix |A]), but now with non-commutative 
skeleton weights iyi(m). 

In particular, the non-commutative "redundant" weights Wjj-, i— j > 1, generalizing (|A.ip now 
read, as a consequence of iterations of ()3.5p : 



(3.7) 



= i ( n y2ay2a-l ) 
\a=i-l J 



where the product is taken in decreasing order, a = i — l,i — 2, + 1. 

Example 3.3. In the case of the fundamental Motzkin path iuq = (0,0, ...,0), we have 

(3.8) F^,{t) = l + a(")(tyi,ty2,...,ty2r+i)yi 

with Yi = yj(mo). The canonical form is obtained by applying (j3.6p at level zero, with a = 0, 

b = 1 and c = tyi, with the result: 

(3.9) 

Fmo(i) = fl - i(l - t{l - tys - t(- • • (1 - ty2r-l - t{l - ty2r+l)"V2r)"^ • • • )"V4)"V2)"Vl 



This is the generating function for paths on the graph of Example \A.^ from vertex 1' to 

itself, with the non- commutative weights Wj^j/ = Wj_i^j = 1, Wj/^j = ty2i-i, = ty2i-2 for 

i >2 and wi'^i = 1, wi^i/ = tyi. 

Example 3.4. In the case m = mi = (0,l,2,...,r — 1), we have mi = rrii-i + 1 for all i. The 
corresponding continued fraction reads, with yi = yi(mi).- 

(3.10) F^,it) = 

1 + t (l - tyi - ty2 - t^(..(l - ty2r-l - ty2r - t^(l - ty2r+l)~V2r+iy2r)~"^..) V2yi) yi 



^1 - t{l - t(..(l - t{l - ty2r-+l)"V2r-)"'..)"V2) Vl 

where the first line is the Jacobi-like expression (j3.4p and the second line has the canonical form. 

The second identity in Example 13.41 is the natural non-commutative generalization of the Stielt- 
jes continued fraction expansion. In the following, we will refer to the case m = mi as the "Stieltjes 
point" . 

The graph associated to this example is a chain of 2r -|- 2 vertices labeled 1, 2, 2r -|- 2 
(see Appendix E]) . with weights Wj^j+i = 1 and vi^i+i^j = tyi, i = 1,2, 2r -|- 1. 

Example 3.5. In the case r = 3 and m = (0, 1,0), we have, with yi = yi(0, 1,0).- 

F(o,i,o)(*) = 

1 + 1 (l - t(yi +y2)-t^{l- tya + y^ V4 - (1 - tys - t(l - ty?)" Ve)" V5" V4) " V3y2) yi 

\ _ t{l - t{l - t{l - t{l - tys - t(l - tyrr^erHy^ + (1 - ty7)"V6y5"V4))"V3)"V2)"Vi 



-1 
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where the first expression is the Jacobi-like form (|3.4p . and the second is the canonical form, 
obtained by first applying (jS.Sp with a = tya, b = ty^ and c = ty^ while u = 1 — ty^ — t{l — 
ty7)~^ye, then (j3.6p with a = tyi, b = ty2, c = ty3, and finally (j3.6p at "level zero", with 
a = 0, b = 1 and c = tyi. We see the emergence of the "redundant" weight ty^y'^^y^^. See the 
corresponding graph r^o^i^o) AppendixlAi 

3.2. Mutations. In the non-commutative case, we do not have an evolution equation for cluster 
variables. Instead, we start by defining mutations of the weights y. We will then show that 
these mutations have the same effect on the generating functions Fm{t), considered as continued 
fractions, as the non-commutative analogues of the rearrangement Lemmas 12.121 and 12.131 

We consider the set of Motzkin paths m S M. We use the same notion of mutations on 
Motzkin paths as in the commutative case. Recall that mutations act on Motzkin paths as 
follows: /Uj(m) = m' where m' differs from m only in that m[ = rui + 1. In order to obtain 
any Motzkin path m E M via a sequence of mutations, we need only consider the two types of 
mutations described in cases (i) and (ii) of Section 12. 4[ 



(3.12) Case (ii): < 



3.2.1. Mutations of weights. We now define the action of mutations on A, as follows. 
Definition 3.6. Let ^j(m) = m'. Define y' := y(m') as a function of y := y(m) as follows: 

y'2i-i = y2i-i + y2i, 

(3.11) Case (i): { ya^ = y2i+iy2i(y2i-i + y2^)"^ 

y2i+i = y2j+iy2i-i(y2i-i + y2i)""^- 

y2i-i = y2j-i + y2i, 

= y2j+iy2i(y2j-i + y2i)^"'^, 

y2i+l ~ y2i+iy2i-l(y2i-l +y2j)~"'", 
y2j+2 — y2i+2y2i-l(y2i-l +y2i)~"'"- 

with all other y'- = yj unchanged. 

Our main goal in this section is the proof of Theorem 13. 101 concerning the action of mutations 
on the generating functions Fin(t). Towards this end, we prove some non-commutative analogues 
of the continued fraction rearrangement lemmas of Section [2j 

3.2.2. Mutations as non-commutative rearrangements. In addition to the two identities ()3.5p and 
()3.6p . we have the two following local rearrangement lemmas for non-commutative continued 
fractions. 

Lemma 3.7. For all a, b G yi we have 

(3.13) 1 + (1 - a - b)-ia = (l - (1 - b)-^a)"^ 

Proof. By direct calculation. This is the "level zero" application of identity p.6p . □ 
Lemma 3.8. For all a, b, c, d € A, with a + b invertible, we have 

(3.14) a + (1 - (1 - d)-ic)"^b = (1 - d - (1 - c')"^b')"^ a' 
where a',b',c' are defined by 

(3.15) a' = a + b b' = cb(a + b)^ c' = ca(a + b)"^ 

Proof. By direct calculation. □ 



Together with the identities (|3.5p - (|3.6p . these two lemmas may be used to rearrange the non- 
commutative continued fraction expression of Fin(t) ^ ^^/{t) (|3.4p . 

Define the non-commutative version of the fractions Jj as in the proof of Lemma 12.131 
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Definition 3.9. The non- commutative finite continued fraction Jj(ni) is defined inductively as 
follows: 

Ji(m) = y2i-i + (1 - Ji+i(m))"^y2i, Jr+i(m) = y2r+i(m). 

Note that Jo(ni) = J''''^(y(m)), and Fni(i) = 1 + tJo(m)yi, with the convention that y_i = 
and yo = 1- 

Let us fix i, and define: 

a = ty2i_3,b = ty2i-2,c = ty2i-i,d = ty2j,e = ty2i+i,f = ty2i+2,u = (1 - Jj+2(m))"\ 

where y^ := yj(m). Similarly, we use the same notation for the primed counterparts y^ := yj(m'). 

Starting from Fm(t) in Jacob i-like form (j3.4p . we apply fraction rearrangements on the part 
J,_i of F^{t). 

• Case (i): We first transform the Jacobi-like form into the local canonical form for Jj_i(m) 
by applying ()3.6p at the second level, namely: 

Jj_i(m) = a+ (l -c-d- (1 -e-uf)"^ed)"^b 
= a+ (1 -c- (1 - (1 -uf)-ie)-M)"^b 
Next, we apply Lemma 13.81 to the second level, namely: 

Ji_i(m) = a + (1 - (1 - (1 - e' - uf)-M')"^c')"^b 

with c' = c + d, d' = ed(c')^^ and e' = ec(c')^^. Finally, we put this back in Jacobi-like 
form, by applying (j3.6p backwards to the first level, namely: 

Ji_i(m') = J»_i(m) = a + b+ (l - c' - (1 -e' - uf)"M')"^c'b 




• Case (ii): We apply directly Lemma 13.81 to the second level of Jj_i(m), namely: 

Ji_i(m) = a+ (1-c- (l-e-uf)"M)"^b 

= a + (1 - (1 - (1 - e' - uf')"^(d' + uf'(e')"M'))"^c')"^b 

with c' = c + d, d' = ed(c')^-'^, e' = ec(c')^-'^ and f = fc(c')~^. To bring this back to 
Jacobi-like form, we first apply ()3.5p backward on the third level, and then ()3.6p backward 
on the first one: 

J,_i(m) = a+(l-(l + (e')-M'-(l-e'-uf')-^(e')-M')-^c')"'b 
= a + b+(l-c' + (e'y^d' - (1 - e' - uf')~^(e')"M')"^c'b 
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Note that in both cases, the relation of the primed variables to the unprimed variables is the same 
as the relation of the mutated variables y'^ to yj in Definition 13.61 

Clearly, for i > 1, we always have Jo(m') = Jo(m), or equivalently Fm(i) = FmKO- When 
i = 1, we have a = 0, b = 1 and in all cases c' = yi(m'), so the final formula amounts to 
Jo(m) = H-tJo(m')yi(m'), or equivalently Fm(t) = I + tF^,{t)yi{m). 

We therefore have the following theorem: 

Theorem 3.10. Under the mutation fii on A, fii : y(m) i-^ y(m') of Definition \3.(^ the function 
Fm(i) is transformed as: 

(3.16) i = l: ¥^{t) = l + tF„,(t)yi(m), 

(3.17) i>l: ¥^{t) = ¥^,{t). 

3.3. An explicit expression for the weights. Although we do not have evolution equations of 
the form ()2.ip in the non-commutative case, we make a definition of the variables R„ as follows. 

Definition 3.11. Let RnR^"^ he defined as the coefficients in the power series expansion of (t), 
as follows: 

oo 

(3.18) F™o(i) = ^rR„Ro\ 

n=0 

In particular, yi(mo) = RiRq^ . To make the contact with the commutative case of Section 2, 
one should compare R„ with of that section. 

In this section, we relate the weights y(m) to the coefficients of the generating series The main 
result is Theorem 13. 331 below, which gives an explicit expression for all y(m) in terms of the Rn's. 

3.3.1. Stieltjes point and quasi- determinants. The generating series Fmo(i) may be written as the 
non- commutative "Stieltjes" continued fraction Fjn^(t) in its canonical form (j3.10p at the Stieltjes 
point m = mi. Since the sequence of forward mutations from mo to mi only involves mutations 
m with i > 1, the generating series is preserved at each step, and 

Fmo — FjT^j . 

The coefficients yi(mi) of the non-commutative Stieltjes continued fraction (j3.10p are expressed 
in terms of the Rn's via a non-commutative generalization of the celebrated Stieltjes formula 
(Theorem 13.171 below). For this, we need a number of definitions. In the following, A denotes a 
square n x n matrix with coefficients ajj € A. 

We recall the notion of quasi-determinant [18] . 
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Definition 3.12 (Quasi-determinant). The quasi- determinant of A with respect to the entry (p, q) 
is defined inductively by the formula: 



a, 



q\. .] 



Oil 



where A^'^ denotes the n — 1 x n — 1 matrix obtained from A by erasing row p and column q, while 
keeping the row and column labeling of A. 



We refer to |18j for a review of the many properties of the quasi-determinant. For now we will 
only need the following: 

Proposition 3.13 ([18J, Proposition 1.4.6). The three following statements are equivalent: 

• (i) \A\p^q = 

• (ii) the p-th row of A is a left linear combination of the other rows 

• (Hi) the q-th column of A is a right linear combination of the other columns 

Note that the quasi-determinant is not a generalization of the determinant, but rather of the 
ratio of the determinant by a minor: 

Proposition 3.14. If A is commutative, and the minor \AP'''\ ^ 0, then the quasi- determinant 
\A\p.q reduces to a ratio of the determinant of A by its p,q minor: 



(3.19) 

Finally, we also need: 



1^ 



1^1 



\AP'i\ 



Definition 3.15 (Quasi- Wronskian). The quasi- Wronskians of the sequence (RA/)MeZ o'^e defined 
for i > 1 and n £ Z as the following i x i quasi- determinants: 



(3.20) 



A,; 



R-n+i-2 R-n+i-S 



H-n-l 



Rn- 



1,1 



Example 3.16. We have for i = 1,2, 3: 



Ai,„ 


= Rn 


A2,n 


= Rn 


A3,n 


= Rn 



Rn 



-R„fRn-l — RnR^ix-'^"-2) ^Rn+1 — Rn+l(Rn — Rn-lR„i2-'^"-l) """^^n+l 



We are now ready for the non-commutative Stieltjes formula [19]. We start from the (infinite) 
non-commutative continued fraction S(t;x) with coefficients Xi € A, i € Z>o- It is defined 
iteratively by 



(3.21) 



S(t;xi) 
S(i;xi,X2, ...,Xfc) 
S(t;x) 



(l-txi)-i 

S(t;xi, ...,Xfc_2, (1 - txfc)~^Xfc- 
lim S(t;xi,X2, ...,Xfc) 
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Theorem 3.17. The coefficients Xj of the non- commutative continued fraction S{t;x) are related 
to the coefficients of the corresponding power series S(t;x) = ^„>q t"R„RQ ""^ via: 

(3.22) X2^_i = Ai^iArli 

(3.23) X2i = Ai+i,iAr/ 
for all i > 0. 

Proof See Ref. [19J, Corollary 8.3. □ 

Comparing S(t;x) (|3.21|) with the finite canonical continued fraction at the Stieltjes point 
Fmi (t) ()3.10p , we see that the latter is a finite truncation of the former, which may be implemented 
by imposing that X2r+2 = 0, while Xj = yi(mi) for i = 1,2, 2r + 1. From Theorem 13.171 this 
finite truncation condition is equivalent to the equation Ar.+2,T-+i = 0, and Theorem 13. 171 gives us 
an evaluation of the non-commutative weights at the Stieltjes point in terms of the Rn's. In fact, 
we have: 

Theorem 3.18. The coefficients R„ of the power series expansion of the finite continued fraction 
Fmi(i) satisfy: 

(3.24) A,+2,n = (n>r) 

Proof. As a consequence of the definition (j3.1Up . we may write (with y-j = yj(mi)): 

Fmi (t) = P{t; yi, y2r+i)~'Q(t; y2, y2r+i) = P{t)-^Q{t) 

with polynomials P,Q G A[t], such that P{t),Q{t) = l + 0{t) and deg(P) = r + 1 and deg(Q) = r. 
Indeed, from the inductive definition (j3.2ip we get 

(3.25) P{t; yi, y^+i) = P{t; yi, ys) - tys+iP{t; yi, y^-i) 

while P(t;yi) = 1 - tyi, and Q(t; y2, y^) = P(t; 0, y2, y^). Writing P(i; yi, y2r+i) = 
'^l=oi—tyPi, we conclude that the Rn's satisfy a linear recursion relation of the form 

r+l 

(3.26) 5^(-l)'P,R„+,+i_, = 

i=0 

for all n > 0, with Pq = 1 and Pr+i = y2r+iy2r-i • • 'yi- This implies that the first row of the 
Wronskian matrix of (j3.20p is a left linear combination of its other rows, and the theorem follows 
from Proposition 13.13] □ 

Remark 3.19 (Conserved quantities). Note that the coefficients Pm, rn = 0,1,..., r + 1 o/ the 

linear recursion relation ()3.26p are conserved quantities under the evolution of the weights of 
Definition \3.(j\ Indeed, throughout this evolution, the denominator of the generating function 
Fm = Dm^m (where both factors are equal to 1 modulo t) remains invariant. Therefore, = 
P{t;yi, ...,y2r+i) defined above, is the generating function for the r+l conserved quantities Pm 
(with Pq = 1). Note, however, that in general these quantities do not commute with each other. 

Remark 3.20. The conserved quantities Pm have an interpretation in terms of non- commutative 
hard particle partition function on a chain, as in [7] in the commutative case. That is, Pm is 
the sum over configurations of m hard particles on the integer segment [l,2r + 1], with a weight 
UimVim-i ' ' ' Vh P^'^ particle configuration occupying vertices ii,i2, ■■■,im such that ij+i — ij > 1. 
This interpretation is a consequence of the recursion relation (|3.25p . 

More generally, one can show that for any Motzkin path m, the conserved quantities, when 
expressed in terms of the weights yj(m), are (non-commutative) hard-particle partition functions 
on certain graphs introduced in 0. 
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For later use, it will be convenient to extend the definition of R„ to n < 0. This is easily done 
by extending the linear recursion relation (13.26P to any n € Z. For instance, this allows to define 
R_i = P~_^i ^l^oi—^Y'^^ Pi^r-i- The main equation extends accordingly to 

(3.27) Ar+2,n = for aU n € Z 

Remark 3.21. If A is commutative, and Ri^n o-i"s solutions to the Ar Q-system (|2.ip . we have 
R„ = Ri^n = Rn- Then according to Proposition \3. 14\ we have: 

^ _ det {{Rn+i+l-a-b)l<a,b<i) _ Ri,n 

det {(Rn+i+l-a-b)2<a,b<'i) Ri-l,n-l 

The finite truncation condition (j3.27p is equivalent to i?r+2,n = 0, for all n (z 

3.3.2. Properties: from Pliicker to Hirota. In this section, we derive a recursion relation for the 
quasi- Wronskians, which is a non-commutative analogue of the discrete Hirota equation. To this 
end, we must first review a few useful properties of quasi-determinants [18]. 

We recall the so-called heredity property of quasi-determinant (see [T8|, Theorem 1.4.3). This 
states that we get the same result if we apply the definition 13.121 to any r x r decomposition of 
the matrix A into possibly rectangular blocks, viewing the blocks as some new non-commutative 
matrix elements, and then evaluate the quasi-determinant of the resulting matrix. This holds 
only if the {p, q) entry of A belongs to a square block, and all expressions (involving inverses) are 
well defined. 

Using this property, we immediately get the following two lemmas: 

Lemma 3.22. Let A be a k+lxk+1 matrix with last column such that ai^k+i = fori = 1, 2, k. 
Then for 1 < p,q < k, we have: 

(3.28) \A\,,, = \A'+'''+X,, 

Proof. By the heredity property, we decompose the matrix A into a first k x k square block 
B = ^^+1.^+1 which contains the {p, q) entry of A, a zero k x 1 column vector z, a 1 x k 

B z 

row u = (afc+1,1, afc+i,fc) and the single matrix element x = Ok+i^k+i, to get: = 

\B — ZX~^u\p^g = \B\p^q. □ 

Lemma 3.23. Let A be a k + 1 x k+1 matrix with first column such that Cj^i = for i = 1,2, k. 
Then we have: 

(3.29) l^lfc+1,1 = (i-k+i,! 

Proof. By the heredity property, we decompose the matrix A into a zero k x 1 column z, the 
kx k block B = yl'^+^'i, the single matrix element x = a^+i i, and the 1 x k row vector u = 
(afc+1,2; •••5 afc+i,fc+i): 



\A\ 



fc+1,1 



z B 
X u 



= \x -uB ^z\k+i,i = X = afc+1,1 

fe+1,1 

□ 



Quasi-determinants and minors are related via the so-called homological relations: 
Proposition 3.24 (Homological Relations [18], Theorem 1.4.2). 

(3.30) ({A'^'keY' \Akj = - \A\k,j 
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Quasi-determinants obey the following non-commutative generalizations of the Pliicker rela- 
tions. Let A he a k X m matrix with elements in A, and m > k. Let A{ii, ...,in) denote the k x n 
matrix made of the columns of A indexed zi, z„. 

Definition 3.25. For any ordered sequence I of column indices of A, we define the left quasi- 
Pliicker coordinates 



(3.31) 



mi, I) 



,s,ir'\A{j,i)\sj 



Note that the definition is independent of the choice of row index s. Moreover, if j is among /, 
then q- j = 0. 

We may now write the generalized Pliicker relations: 

Proposition 3.26 (Non-commutative Pliicker Relations, [18J, Theorem 4.4.2). Let M = mi,m2, mfc_i 
and L = ii,i2, ■■■,^k be two (ordered) sequences of column indices of A, then for any i L: 

3&L\M 



(3.32) 



With these properties in mind, we now turn to our main theorem: 

Theorem 3.27 (Discrete Non-commutative Hirota equation). The quasi- Wronskians obey the 
following relations 

(3.33) ^i+l,n = ^i,n+l — '^i^ni^i^n-i — A._\^„)Aj^„ 

for i > 1 and n E Z, with the convention Aq ^ = 0. 

Proof. Note that using the initial values i = 1,2 of Example 13.161 we get A2,n = R-n+i — 
Ita^n-i^n ~ ^i,n+i ~ ^i./t^Tn-i^i,"' which agrees with (I3.33P for i = 1, given the con- 
vention for Aq^ = 0. We now turn to the general proof of the identity (|3.33p for i > 1. Let us 
apply the Pliicker relation (I3.32p to the following i + 1 x i -|- 3 matrix and column index sets: 



/ 1 R-n+i R-n+j-1 
R„+i-l R„+.j_2 



A 



Rn 

Rn-1 







Rn+1 
Rn 



R« 



M = (3,4,...,i + l,i + 2) 
L = (2,3,...,i + l,i + 3) 



Rn-i+l 
\0 R„ Rn-l ••• Rn-i 1/ 

so that L\M = 2,i-|-3. We now compute the various Pliicker coordinates involved in the sum 
(j3.32p . Introducing 
(3.34) 



we find 



R-n+1 



Rn+l 



R-n-1 



Rn-i+2 
R-n-i+1 1 



1,1 



R„+i-i R„+i_2 

R„+i-2 Rn+i-3 



Rn 

Rra-1 



Rn-1 
Rra-2 



Rn 

Rn-1 



Rn-i+1 
R-n— ?' 



1,1 



,A{2} 



9l,2 — ^i+l,n 92,1 

SO the Pliicker relation turns into 
(3.35) 



j-i 



9l,j+3 — '^i+l,n 



L\{i+3} _ -1 



+1 
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By direct application of Lemma [3221 we compute = Aj „+i. To compute '^i+i^n we ap- 

ply the homological relation ()3.30p with i = j = \, k = i + \, li = 2 and with the matrix 
A with entries: a^^^i = and a^^y = B^^-x-y+i for y > I. We get tpi+i^n = 1^1 1,1 = 

-|A*+i'i|i_2(^|Ai'^|i+i,2) The first factor is |^*+^'^|i,2 = ^i,n by definition, while the 

last factor is = Cj+i^i = 1 by applying Lemma 13.231 with k = i. We are left with the task 

of computing lA"*^'^ |j+i^2 = V'i,n-i) where 



1pi,n 



R-n+i-2 
R-n+i-3 



R-n-l 



R 



■n—l 



Rn- 



i,l 



Kn-x-y+i+i, i = i, j = k = l and 



Again, we use the homological relation ()3.30|) with ax,y 
to get: 

V-i.n = -\A^'^\i,2 (|^''^|l,2)~^ l^ll.l = -i^i-l,n~l^l}i^n^i,n 

Recalling that (fi+i^n = — Ai,nV'i"n-i' ^^^^ turns into the recursion relation 

(3-36) ^i+l,n = - (^^i,n^i,n~l) fi,n~l 

Collecting all the results, we re-express (j3.35p as: 

(3-37) ^i+l,n = ^i,n+l — ^i+l,riVi+i^n+l^i,'n.+l 

(3.38) = + <^i+l,n<^i)^\,n 

We deduce that 

A.,„(Ar_\^„ - Ari_^) = A,,„Ar\_ JA,,„_i - A,_i,„)Ar^i_, 



(A 



i+l,n 



H,n+1 



lA; 



and the theorem follows. 



□ 



Remark 3.28. Theorem 3.21 is a discrete analogue of the continuous formula of Theorem 9.7.2 



of [IBJ- Note that the continuous formula contains a misprint (there should he no —1 power for 
the term in parentheses). 

Remark 3.29. In the case when A is commutative, we write R„ = and define for n G Z; 
ro,n = 1, Ti,n = Rn o,nd ri^n = <^Q^i<i,j<i Rn+i+j~i-i- Then, according to Proposition \ 3.14\ we 



have A,- 



j— l,n— 1 ■ 



It is easy to check that in this case the Hirota equation (|3.33p is 



equivalent to the so-called A00/2 Q-system: 

„2 



~l~ Ti+l,n'^i~l,n 



(i > 0,n G Z) 



Further restricting to rr+2,n = for all n & Z amounts to the finite truncation condition ()3.27p . 
We may then pick the boundary condition rr+i,n = 1 for which rj^„ = Ri^n, the solution of the 
Ar Q-system. 



Remark 3.30. The truncation equation Ar+2,n = (I3.27P implies the existence of a conserved 
quantity. Indeed, writing the Hirota equation in the form (j3.37p for i = r + 1, we get that 
(Pr+2,n+i = Vr-+2,n ^-5 independent of n, due to the truncation equation. 
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3.3.3. General expressions for the weights. By applying the non-commutative Hirota equation 
()3.33p . we may compute explicitly all the weights y(m), m G M, as given by Theorem 13. 10^ in 
terms of the Rn's. Let us first define new variables for i >0 and n G Z: 

(3.39) A,n = Aj+i_„A~^ Ci^n = (-l)Vi+l,n 

with the initial conditions -Do,n = and Co,n = 1- Note that (j3.36p implies the relation 

(3.40) Ci^n = Aj^^A- 
hence we have 

i-l 

(3-41) = n A.-,,„-, Ari^. „_ 

j=0 



Note also that Remark 13.301 implies that 

(3.42) Cr+l,n = Cr+1 

is independent of n. 

Remark 3.31. In the commutative case R„ = Rn, with ro,n = 1? ?'i,n = Rn o-nd ri^n solution of 
the ^oo/2 Q -system of Remark \3.29l we have 

i,n — ^i,n — ^i,n — 

'"i— l,n— 1 ^i,n—l ^i,n+l'^i,n 

If we further restrict to the case of the Q -system (by imposing rr+2,n = and rr+i^n = Ij; we 
note that the variables Ci^m and Di^m reduce to the variables Ci^n o,nd di^n of Theorem \2.14\ 

The variables C, D obey the following relations: 

Lemma 3.32. We have 

(3.43) Di^n+l = Ci+l,n+lCi^n'^i,nCi-l,nCi^n+l 

(3.44) Ci,„+iCj_\_„ = A,n + C'i,nC',;-l,„ 

(3.45) Cj+i_„C-_^_^ = Din + C'i+l,nC'i,ri 

Proof. Eq. p.43p holds by definition, as a consequence of (j3.40p . Eq. (j3.44p is equivalent to the 
Hirota equation in the form (j3.38p . upon multiplication by Using (j3.44p and (|3.40p we 

compute 

Di^n -\- Ci+i^nC^ n = Aj+i^nAj_^\^„_]^(A,?i-l + Ci,n-lC'i_\,„_i)Aj,„_lA-^^ 

= Aj+i^„A.^\^,^_-,^Ci,„C-_\^„_;^Aj_„_iA.^^ = Ci+i^nCin_^ 
and the lemma follows. □ 

We are now ready for the main theorem of this section: 
Theorem 3.33. The weights {yj(m)} are expressed in terms of the R„ 's as: 

(3.46) y2i-i(m) = Ci,ni,+iC-\^^_^^^ 

y2i(m) = if = mi+i + 1 

* 1 otherwise 

(3.47) X X 



Ci-i,m,+iCj_\,m^_j+i if m - mi_i 



1 otherwise 
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Proof. As already pointed out, the recursions (|3.1ip - ()3.12p determine uniquely the weights y(m) 
in terms of the initial data y(mo) at the fundamental point mo = (0, 0, 0). However, Theorem 
13.171 with X2r+2 = gives expressions for the weights Xj = yj(mi) at the Stieltjes point mi = 
(0, 1,2, ...,r — 1). So we will assume that eqns. (|3.46p - (j3.47p hold at the fundamental point, and 
prove by induction that they hold everywhere. We then simply have to check that ()3.46p - ()3.47p 
hold at the Stieltjes point and the theorem will follow. 
Using (106]) - ([071) for rui = i - I, we get: 

y2j-i(mi) = Ci^iC~\^^_^ y2j(mi) = A,j 

These reduce respectively to by (13:1(1 and to (13:231) by (lOOD . So the solution ([336])- ([337])) 

holds at the Stieltjes point. 

We are left with proving the inductive step: assume (I3.46p - (l3.47p hold at some point m, with 
weights Yi = yj(m), then show that it does at m' = /x^m, with weights y- = yj(m'). We must 
distinguish the two usual cases: 

• Case (i): mj_i = rrii = m, rrij+i = m + 1. By (I3.46p - (|3.47p we have: 

y2j-l = Ci,m+lC'i-l,m+l' = Di^m+li Y2i+l = Cj+l,m+2C'j^m+l 

The evolutions ()3.1ip give: 

y2i-i = y2i-i + y2i = Ci^rn+2C~\^j^i (by eq. (|3.44p) 

y2i = y2i+iy2i(y2i-i)~"^ = A,m+2 (by eq. p.43p) 

y2i+l = y2i+iy2i(y2i-l) ^ = C'i+l,m+2Cj,m+2 

and (j3.46p - (|3.47p hold for m', with m^_i = m and m'- = m'-^^ = m + 1. 

• Case (ii): mj_i = rrii = mj+i = m. Now we have by (l3.46[ )- (13.47p : 

y2i-l = Ci^rrt+lC'^_\^m_|_i> y2j = -Dj^m+l, y2j+l = Ci+l,m+lCj.m+l' y2j+2 = ^A+l,m+l 

where A = Cj+i^m+iCTm+i "^i+2 = m — 1, and A = 1 otherwise. The evolutions (I3.12P give: 

y2i-i = y2i-i + y2i = Ci,m+2Cil\,„+i (by eq. ^M) 

y'2i = y2i+iy2i(y2i-i)"^ = C'i,m+iC^\,^+2A,m+2 (by eq. (|3.43P) 

y2i+l = y2j+iy2i-l(y2i-l) ^ = C'?;+l,m+lCj_m+2 

y2i+2 = y2j+2y2i-l(y2j-l) ^ = ^A+l,-m+lCj,m+lCj^Jj^_2 

and (j3.46p - (|3.47p hold for m', with rn[_i = m'^^^ = m and m'^ = m + 1. □ 

Example 3.34. The particular set of non- decreasing Motzkin paths with mi < m2 < • • • < m-r 
is obtained by only applying mutations pertaining to the case (i). The corresponding weights are 
particularly simple: 

y2j-l(m) = Ci_mi + lCjl\,m._i+i, y2i(m) = A,mi + 1- 

In particular, for the fundamental point mg, we get 

y2j_i(mo) = Ci,iCjl\ 1, y2j(mo) = 
This can he compared to the so-called J -fraction expressions of Theorem. 8.12. 
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Example 3.35. In the case r = 2, we have the three Motzkin paths (0,0), (0, 1) and (1,0). We 

have Co,n = 1 and C-^^n = C'a- The weights are: 



yi(o,o) = 




yi(o,i) 


= 


yii 


:i,o) 


= Ci^2 


y2(o,o) = 




y2(o,i) 


= 


y2i 


:i,o) 


= C2,lC2 2-^1^2 


y3(o,o) = 




y3(o,i) 


= (^2,2(^1 1 


yai 


:i,o) 


= C'2,lC^2 


y4(o,o) = 


D2,l' 


y4(o,i) 


= I?2,2 ' 


y4i 


:i,o) 


= L'2,iCi^iC]^ 2 


y5(o,o) = 




y5(o,i) 


= (^3(722 


ysi 


:i,o) 


= CZC21 



In the formulation as paths on Tm via the canonical form of the continued fractions ^ the only ^^re- 
dundant" weight occurs in the last case m = (1,0), where W3^i = ty4y3^"'^y2 = i-D2,iCi_iC^2-^i^2- 
Note that weights may involve R„ 's with n < in their definition (e.g. Aj^„ for n < i — 1). 

Theorem 3.36. We have 

ra>0 

In other words, the quantity Yii^n+mi^i\nx partition function of paths from vertex 1 to itself 

on Tin, and with the non- commutative weights (j3.ip associated to (|3.46p - (|3.47p . ^45 such it is a 
Laurent polynomial of these weights, with non-negative integer coefficients. 

Proof. The statement is proved by induction on mutations. It holds for m = mo, with mi = 
(mo)i = 0. Any given m G M is obtained by an iterative application of mutations fii of type (i) 
or (ii). Assume the theorem holds for some m € M. Let us form m' = ^^m. When i > 1, we 
have Fmit) = Fm'it) = Z]n>o *"^"+miR'ml by (|3.17p : hence the same quantity Kn+mi'R-m\ is 
interpreted as the partition function for paths on T^/ with the weights associated to yj(m'), and 
the theorem holds for m', as m'l = mi. If i = 1, then we have p.l6p . Using yi(m) = Ci^mi+i = 
R-mi+iR-ml, we get 

Fm(i) = ^ i^R-n+miH-ml = 1 + tF^/{t) Rmi+lR-ml 

ra>0 

Fm'(*) = ^""R-n+mi+lR-m'l+l 
ra>0 

and the theorem holds for m', with m'l = mi + 1. □ 

Example 3.37. Let us apply Theorem \ 3. 36\ to the calculation of R3 in the case r = 2, at the 
points m = (0, 0), (0, 1), (1,0) respectively. We have 

m = (0, 0) : R3R0 ^ = (yi + y2) Vi + (ys + y4)y2yi 

= (Ci,i + ^i,i)'Ci,i + (C2,i + D2,iCi,i)C^lDi^iCi,i 

m = (0, 1) : R3R0 ^ = (yi + y2)^yi + y3y2yi 

= (C'l,! + ^l,l)^Cl,l + C2,2CilDi^iCi^i 

m = (i,o): R3Rr^ = yf + y2yi + ys.iyi 

= (^"1,2)^ + {C2,l + D2,lCl^l)C2^Di^2Cl^2 

respectively equal to the partition functions for non- commutative weighted paths from the root to 
itself on the corresponding graphs Fm, with weights = yj(m) given in Example \3.35\ and with 
respectively 3, 3, 2 steps towards the root. The identities between these expressions are obtained 
from identities of Lemma \3. 3S\ relating C and D: 

Cl,l + Dl,l = Ci^2, C2,l + -C>2,lCi_i = 6*2,2, C^^Di^i = 6*2 2-01,2^1,2 
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4. Application I: T-system 

Here, we apply the construction of the previous section to the Ar T-system (ll.ip . which can be 
written as a non-commutative Q-system as shown in [9]. The boundary condition implies a 
truncation relation Tj.^2,j,k = for all j, A; G Z, similar to that studied in Section [3l so we can use 
the results of Section [3] to give an alternative derivation of the T-system positivity results of [9]. 

4.1. The Ar T-system as non-commutative Q-system. We may view the T-system equations 
(jl.ip as matrix elements of a relation in a non-commutative algebra in a particular representation. 
We consider the non-commuting variables Tj ^. € A which have a faithful representation on a 
Hilbert space "K with basis {\i)}i^z- 

(4.1) T,^k\j + k + i) = T,,j,k\j -k-i) {i>0;j,keZ) 

Introducing the shift operator d which acts on the same representation as d \t) = \t — 1) we may 
rewrite the Ar boundary conditions Toj^fc = T^+ij-,fc = 1 as: 

(4.2) To,fc = d2^ T,+i,fe = d2('=+'^+^), 
while the finite truncation condition implies 

Tr+2,k = {kez) 

Note that all variables Tj fc, < i < r -|- 1, are invertible for generic non- vanishing matrix 
elements Tjj- fc, with 

(4.3) iT^,kr^\j-k-i) = -^\j + k + ^) {0<i<r + l; j,keZ) 
We have the following 

Proposition 4.1. The Ar T-system (jl.ip is obtained as matrix elements in "K of the following 
non- commutative Q-system: 

(4.4) Ti,fe+i(Ti,fe)-iTi,fc„i = Ti^k + T,+i,fc(T,,fc)-iT,_i,fc (1 < i < r; k e Z) 
with the boundary conditions ()4.2p . 

Proof. Applying (14. 4p to the vector \j + k — 1 + i) and using (14. ip and ()4.3p yields (II. ip divided 
by Tjj+i,fc times li - /c 1 - i). □ 

Remark 4.2. Note that the operators Tj^^ are only "mildly" non- commuting, in the sense that 
the operators {d^Tj fcd~'"~^*~^'^'}m,fcGZ o,'^^ M diagonal when acting on Ji, and therefore com- 
mute with each other. In a sense the d operator is the only source of non- commutation. It was 
interpreted in [9] as a time-shift operator, and the non- commutation of the variables is due to 
their time- ordering. Nevertheless, this gives us a first non-trivial example of a non- commutative 
Q-system. 

4.2. Quasi- Wronskian formulation. The Ar T-system is a particular case of the non-commutative 
finite truncation equation (|3.27p . We start from the solution Tj j ^ of the system (jl.ip and the 
operators Tj^ (j4.ip - (|4.3p . This gives rise to the sequence R„ = Ti^„ G ^l, n € Z. Defining 
^i,n, n E Z as the quasi- Wronskians (j3.20p of this sequence, we have the following action of the 
quasi- Wronskians on the basis {\t)}t£z- 

Theorem 4.3. 

(4.5) A,,„ \j + n + i)= \j-n- i) (1 < i < r + 2; j,n e Z) 
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Proof. By induction on i. For i = 1, we have Ai^„ = R„, and Rnb + n + 1) = Tij^^ \j — n — 1), 
which amounts to (I4.5j) for i = 1, as ?oj-^fc = 1 for ah j,k. Assume (I4.5P holds for all /3 < i. 
Using the Hirota equation ()3.33p . and the recursion hypotheses we get: Aj+i „|j + n + i + 1) = 
Si+i,n \j -n-i-1), with 

rj~i / rjl rj-i \ rj-i 

^i,j,n+l J-i,j — l,n I -^i—2,j,n—l -^i— lj',n— 2 \ -^ij+l,?! 



Oi+l,n — 7^ H 



i—l,j,n ^i—l,j,n—l \ -^i— lj,n -^i,j^n—l J J-i—l,j+l,n—l 
Ti—lj^n Ti—ijnT{jn—l 
-^l!,jr',n— 1 

where the second line follows from (jl.ip for z ^ z — 1 and k = n — 1, while the third follows from 

([TT]) with A; = n. □ 

Remark 4.4 (Relation to the known determinant relation of [2]). The T -system (ll.ip may 5e 

simplified by first eliminating the Tij^n, i > ^ and j,n ^"L in terms of the Tj^n = ^i,j,n by use of 
the Desnanot-Jacobi formula, as the following determinant: 



TiJ^n — det(0jj^„), ( 0i,jf,n ) — Tj_aJ^b,n+i+l-a-b (1 ^ 0,,b < i) 



We note that 



_ \r\ I 
— l'-'ij,n|l,l 



2^— 1 j,n— 1 



SO the diagonal operator d " *Aj_„d " * has the eigenvalue |0ij,n|i,i when acting on for all 
j€Z. 

Applying Theorem 14.31 for i = r + 2, we find: 
Corollary 4.5. The T-system solution (Ra;)^^^ satisfies the finite truncation equation (|3.27p . 

We may therefore apply the results of Section [31 The C, D operators (|3.39p act on the basis 

{\t)}tez as: 

a,n\t) = \t-2i) Cr^\t) = ^4±I^±i±l^\t + 2i) 



Ti-\^l^t+n+i—l,nTi—l^t+n+i,n—l 
Ti^t+n+i—l,n—lTi^t+n+i,n 



Using the notation 

/ . „\ Ti^t,m r ^i+l,t,m^j— 1,4+1, m—1 

(,4.bj Ci^t,m — Tf, 0'i,t,m — Tf, 7f^ 

and the fact that expressions of the form Ci^rnCj^p act diagonally, Theorem 13.331 yields non- 
commutative weights yj(m) which act on the basis as: 

(4.7) y2i-i(m)|t) = ^i,t+rn,±izlfIh±^\t-2) 

y2i(m)|t) = <^ c,+i,t+™^^+,,™^ + i ^ \ X di^t+m,+i,m,+l 

I 1 otherwise I 

; ■ II rui = rrii-i — 1 , 

1 otherwise ( 
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Remark 4.6. To he able to compare the results ()4.7p - (|4.8p to those o/[9], we note that here, in the 
non- commutative weighted path interpretation, the up steps have weights 1, acting as the identity 
on the basis {\t)}teZ) while all the down weights (including the redundant ones) act as a diagonal 
operator times (P. Interpreting t as a discrete time variable, we see that up steps use no time, 
while down steps use two units of time. In Ref. the time- dependent weighted path formulation 
of the T -system solutions uses skeleton down steps as well as up steps that all use one unit of time 
(i.e. are the product of a diagonal operator by d), while the redundant steps may go hack in time 
(a step i ^ i — k, k > 1 uses 2 — k units of time, hence the corresponding weight is a diagonal 
operator times d^~^J. The two descriptions are equivalent up to a "gauge" transformation, namely 
a conjugation of the weights (j4.7p - (j4.8p .' y2j-i — > d^*^^y2j-id* andy2i — ?• d~*~-^y2jd'. 



5. Application II: quantum Q-system 

5.1. Quantum cluster algebra. The solutions of Q-system for have been shown to form 
particular clusters in a suitable cluster algebra |2T]. On the other hand, cluster algebras are 
known to have quantum analogues, the so-called quantum cluster algebras [3j. In particular, 
there exists a quantum deformation of the Q-system relations. 

The cluster algebra for the Q-system is a rank 2r cluster algebra with trivial coefficients, 
and with initial cluster xq and exchange matrix Bq: 

Xq = (-Rl,0) -^2,0) •••) -Rr,0) -^2,1, -Rr.l) ~ \ C / 

where C is the Cartan matrix for A^-. The subset of (forward) mutations we have considered so 
far is such that each mutation amounts to one application of the Q-system equation, allowing to 
update the cluster via a substitution of the form: Ri^n-i Ri,n+i = (Rf ^ + Ri+i^nRi-i^n) / Tii^n-i 
(mutation /i, if n is odd, ^r+i if n is even). Recall also that the compound mutation /i = 
fJ'2rfJ'2r-i ' ' ' 1^1 acts ou Xq as a global translation of all the n indices by +2 (i.e. Ri^n Ri,n+2 
for all i), preserving Bq. The "half-translation" of all indices by -|-1 is realized by acting only 
with /ieven = /^r/^r-i • • • /^i or //odd = /U2rA^2r-i ' ' ' fJ-r+i in alternancc, which reverse the sign of the 
exchange matrix. 

A quantum cluster algebra with trivial coefficients corresponding to this involves the same 
exchange matrix, subject to the same evolution under mutation, but the cluster variables are now 
non-commuting variables Rj^ri £ ^ whose commutation relations within each cluster are dictated 
by the exchange matrix as follows. We have for any pair (xi,Xj) of variables within the same 
cluster the following g-commutation relation: 



(5.1) 



where A is a suitably normalized 2r x 2r matrix proportional to the inverse of the exchange 
matrix. Note that in a different cluster the commutation relations are different, as they are 
entirely dictated by the exchange matrix, and may involve other cluster variables. 

For the initial cluster, in order for A to have only integer entries, we choose A = (r + 1)Bq^ 


-A 

variables of the fundamental cluster: 



, where A = (r + 1)C~^. We obtain the following g-commutation relations between 



Lemma 5.1. For 1 < i, j < r , we have: 

(5.2) R,-oRj,i = g^'-^Rj,iR7,o 

(5.3) Ri,riR-j> = R-j.nR-j.n {n = 0, 1) 
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where 

(5.4) A,j = Mm(i,j)(r + l-Max(i,i)) (1 < i, j < r) 

Proof. By explicitly inverting the Cartan matrix C, with entries 26ij — D 



Note that the relations (15. 2p and (15. 3p are preserved under the above-mentioned half-translation 
of all n indices by +1, as Bq —Bq, and therefore A — )• — A for /^even and vice versa for ^odd- 
Hence the relations (I5.2I)-(I5.3D extend to 



(5.5) Rj^„Rj^„_i_i — q'^'-'^Hj^n^iTli^n 

(5.6) R^^T^Rj,?! — Rj,nR2,n 

for all n G Z. In particular, the initial sequence R„ = Ri^m n € has the g-commutation 
relations: 

(5.7) R„R„+i = q'^ R„+iR„ 

In the following we consider the same subset of cluster mutations as in the commuting Ar 
Q-system cluster algebra, leading to a subset of cluster nodes indexed by Motzkin paths m G M. 
Each of these clusters Xm is now made of 2r non-commuting variables {'Ri^mi,^i,m^+i}i=i, ordered 
as in the commuting case. As before, these will play the role of possible initial conditions for the 
quantum Q-system below. 

5.2. Quantum Q-system. The particular form of the quantum Q-system we choose is borrowed 
from the T-system case (14. 4p . and compatible up to multiplicative redefinition of the variables 
with the mutations of the quantum cluster algebra. 

Definition 5.2. We define the quantum Q-system for Ar as: 

(5.8) Ri,fc+i(Ri,fc)"^Ri,fc-i = Ri.fc + Ri+i,fc(Ri,fe)~"^Ri-i,fc (i > 1; /c e Z) 
with the boundary conditions 

(5.9) Ro,fc = l Rr+i,fc = l (keZ) 
Note that the latter imply a relation 

(5.10) Rr+2,fc = {k£Z) 

as Ar+i,r+i = 0. We also have R-i,n = as Ao,o = 0. 

Using the (/-commutations ()5.5p - (|5.6p . we may rewrite the system ()5.8p in a simpler form: 

Lemma 5.3. The Ar quantum Q-system is equivalent to: 

(5.11) 5^'>»Ri,fc+iRi,fc_i = (Ri,fc)2 + R,+i,fcRi_i,fc {i>l;k£Z) 
with Ajj as in (15. 4p . and with the same boundary condition (15. 9p . 

Remark 5.4. The relation (j5.1ip is equivalent to the relevant quantum cluster mutations of [3] 
up to a scalar renormalization of the cluster variables. Indeed, defining the variables Xj^„ = 

r+l \ 

q 4 '''.'Rj we find: 

(5.12) Xi,,.+i = g^- (Xi,fc_i)-i(X,,fc)2 + gi(A.,-i+A..+i) (X,,fc_i)-iX,_i,fcX,+i,fc 

in agreement with the definitions of [3\ . Note that boundary conditions remain the same: Xo,n = 
^r+i,n = 1 for all n € Z. 
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5.3. Quasi- Wronskian formulation. The quasi- Wronskians ()3.20p of the variables R„ = Ri,n 
may be expressed as simple monomials in terms of the solutions of the quantum Q-system (15. Sp - 
We have: 

Theorem 5.5. 

(5.13) Ai,„ = <^-^(-i)(2r-+2-0 R. ^R-i^^^_^ (i > 1, n G Z) 

Proof. By induction on i. Define Ti^n = g"^^*"^)^^''^^"*) Ri,nRr-i,n-i' the r.h.s. of (f5T3]l and 
Oi = ^{i — l)(2r + 2 — i). For i = 1, we have Ai.„ = Ri,n = Ti.n. from the boundary condition 
(j5.9p . Assume ()5.13p holds for all 1 < j < i- We use the Hirota equation (j3.33p to express 

(5.14) Aj4-i_n — Tj+i^n = Ti^n+l — Tj+l^n — ^1,7%^ i^n-lO^i,"'-^ ~ '^i-'^,nWi\^n^i,n 

Moreover, defining Bj^^ = (Rj,n)^ + Rj+i,nR-i-i,n - g^''*Ri,n+iRj,n-i (= due to the quantum 
Q-system equation ()5.1ip ). we find that 

rj-l,n - ^i,n-l = q'"" (^(Ri-l,n-l)^ " -^i-l.n-l) R'i-'^2,n-l^i-^l,n-2 

due to the g-commutation relations (|5.5p - (j5.6p and the identity Aj_i^i_i — Aj_2,j-i — Oj = — a^-i. 
Substituting this into (j5.14p . we get 

A — ai+i-p2 "r>— 1 "R~l 

—q 'ri,nrj^„_i(Rj-l,n-l) ^i^2,n-l^i-l,n-2^i-l,n^i,n 
= 

by use of the g-commutation relations (j5.5p - (|5.6p . □ 

Remark 5.6. Theorem \5.5\ allows to interpret Tii^n o,s a "quantum Wronskian", namely a "quan- 
tum determinant" of the Wronskian matrix used in Definition \3.15\ Indeed, we may write 

Ri,n = g^(3r+2-*) Ai_„Ai_i,„_i • • • Ai,„_i+i {i > 1) 

which is a product of principal quasi-minors of the Wronskian matrix, like in the standard defini- 
tions of quantum determinants [TT] (note however that here the quasi-minors do not commute). In 
the cases i = 1,2,3, using the q- commutations and the quantum Q-system, we find the following 
expressions for the quantum Wronskian: 

Rl,n — Rn 

^2,n = ^''Rn+lRn-^l — Rn 

R'3,n = q^'' "^Rn+2RnRn-2 — (Z^^ ^ (Rn+2Rn-l + I^n+lI^n-2) 

+q'' ^{q^^^ + i)R'n+iRnRn-i — Rn 

Applying Theorem 15.51 for i = r + 2, and using the equation (j5.10p . we get: 

Corollary 5.7. The solution of the quantum Q-system (|5.8p - (|5.9p satisfies the finite truncation 
condition Ar+2,n = for all n G Z. 

We may therefore apply the results of Section [3l Assuming Rj^„ is a solution of (|5.8p - (|5.9p . we 
consider the usual subset of clusters x^, also viewed as initial data for the quantum Q-system. 
Recall that the commutation relations between the quantum cluster variables within each cluster 
Xjn are determined by the exchange matrix B^. Using the explicit expressions for the exchange 
matrix B^ of Ref. [7], it is an easy exercise to compute the (^-commutation relations between the 
cluster variables within any given cluster Xm. We find that: 
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Lemma 5.8. 

(5.15) Rj^j^Rj^n-|_p — (f 'Rj^n+p^i,n 

for all n,p,i,j such that the variables belong to the same cluster ^i^. 

Remark 5.9. Lemma [578\ does not cover for instance the case j = i, p = 2, as Ri,n and Rj,n+2 
never belong to the same cluster x^- Using the quantum Q-system and the valid q- commutations, 
we find for instance 

However, we may apply Lemma \5.8\ to compute the commutation of Ai^n and Aj^rn, using Theorem 
\5.5\ when the four variables Rj^^, Ri-i,n-i; Rj.rrn R-j-i.m-i belong to the same cluster x,^;^: 

A. A. _ „(n-m)(Ai,j_i+Aj,,_i)+r(j-i) A . a, 
'-^i,n'-^j,m — y '-^j,m'-^i,n 

Using the ^-commutation relations (I5.15P we easily get the following: 
Lemma 5.10. The quantities Ci^n and Di^n (13.39P are expressed as: 



In particular, the value of the conserved quantity ()3.42p follows from Lemma 15.101 and the 
boundary condition ()5.9p : 

r(r + l) 

Cr+1 = q ^ ■ 
Note that in this case, this conserved quantity is central. 

Theorem 5.11. Let Rj „ be a solution of the quantum Q-system ()5.8p - (|5.9p and R„ = Ri,„. 

The quantity R„+m^R^^ for n > 0, when expressed in terms of initial data Xm, is the partition 
function for non- commutative weighted paths on the graph from the root to itself with n steps 
towards the root and with skeleton weights: 

(5.16) y2.-i(m) = g^-^ Ri,™,+iR-i .R,_i,„,_,R-i 



y2i(m) 



(5.17) 



R-i+l,m,+i+lR-i+i,m,_^jR-i+l,m,R-i+i,m, + l if - + 1 

1 otherwise 

R-j+l,mi+lR'i^mi+l-^j,mi-'^«-l/nii 

R.i-l,.m,+lRr-l,7n,r^«-l,m»-iI^i^\,m,_i+l i^ = '^i-^ " 

1 otherwise 



that are all non-commutative Laurent monomials of the initial data, times some possible integer 
powers of q. As such it is a Laurent polynomial of the initial data with coefficients in 7j^[q,q~^]. 

Proof. We apply Theorem 13.361 with the weights (j3.33p . We then rewrite all the skeleton weights 
using Lemma 15.101 Finally, the redundant weights (j3.7p are also monomials of the initial data, 
with a possible integer power of q in factor. Reordering and rearranging the terms may only 
produce coefficients in Z^[q,q~^]. □ 

Example 5.12. In the case of non- decreasing Motzkin paths (see Example \3.34\ above), we have 
the following expressions for the skeleton weights: 



^Ri,mj + lRj_l^m._j+lRj-: 

y2i(m) = Rj+i,m^+iRj"^^_^^Rr^^Ri_i,m, 



y2i-i(,mj — < I I 
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where we have used the q commutations ()5.15p and the identity Aj_i^j — Aj_i^j_i + i — 1 = to 
simplify the monomials. For the fundamental initial condition xq = Xmo; this gives: 

y2i-i(mo) = Ri_iRj_\ ^Rj qRj-i^o 
y2j(mo) = Ri+i,iR^i Rj QRi_i,o 
For the Stieltjes point Xm^, we have: 

y2i-i(mi) = Rj^jR^_\^^„iRj_i^j_2Rj^/_i 
y2i(mi) = Rj-|.i^jRj^j^Rj^j^_-,^Ri_i^j_i 

Example 5.13. Let us consider the case r = 2 of Example \3. 35\ Setting R„ = Ri^n O'nd 
Sn = R2,n) have the following commutation relations: 

RfiRn+l = (p' Rn+lRrt S^Sji-i-i = Cp' S^+lSn RnSyj = S^Rn 
SnRn+1 = gRn+lSn RriS„+i = Q S^+iR^ {n G Z) 

The A2 quantum Q-system reads: 

Q'^Rn+lRrt-l = (Rn)^ + S„ 
(pSn+lSn-1 = (S„,)^ + R„ 

and the weights at the various Motzkin paths may be expressed as: 

yi(0,0) = RiRo' yi(0,l) = RiRo' yi(l,0) = RaRr' 

y2(0,0) = SiR^iRpi y2(0,l) = SiR^^Ro' y2(l,0) = Sq ^(SO^R-Ir^ ^ 

y3(0,0) = SiR^ISqIRo y3(0,l) = SaRr^R-oSr' y3(l,0) = S1R1S0IR2I 

y4(0,0) = Sr^So-iRo y4(0,l) = S^^S^^Ki y4(l,0) = So^Sr^lRO^Ra ' 



y5(0,0) = Sr^So y5(0,l) = S^^Si y5(l,0) = Sr^Sc 



where we have used (jS.lSp to eliminate the q-dependent prefactors. The redundant weight of 
the last case m = (1,0) is W3^i = ty4y3^y2 = i(1^^^2^^Q^ ■ Note that the conservation law 
y5(m)y3(m)yi(m) = C-^ = q^ is satisfied for the three Motzkin paths. 

Remark 5.14. One can show that the conserved quantities Pm of Remark \3.19\ commute with 
each other in this case. This suggests an analogy with the definition of classical integrability, 
where the integrals of the motion are in involution. In fact, in the Ai quantum Q-system case, 
the Hamiltonian structure and integrability were shown in [13j . The case of the Ai T -system was 
treated in [12]. 

Remark 5.15. The quasi-Wronskians Aj^„ are "bar-invariant" quantities. Let us define an anti- 
automorphism (called usually the bar involution), which we denote by *, by {ab)* = b*a* for all 
a,b € A, 1* = 1 and q* = q^^ , while (Rn)* = Rn for all n G Z. Then, as a trivial consequence 
of the Hirota equation (|3.33p and of the initial condition Ai^n = Rnj we have that 



6. Other non-commutative systems 



In this section, we consider further examples of non-commutative recursions which have the 
positive Laurent property. In the following, all the variables we consider belong to some skew 
field of rational fractions A. 
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6.1. AfRne rank 2 non-commutative systems. The non-commutative Ai Q-system (jl.Sp cor- 
responds to one of the so-called rank two affine cluster algebras, namely with initial exchange 

matrix Bq related to affine Cartan matrices of rank 2. These read Bq = ^ with be = 4, 

6, c G Z>o. While the Ai case corresponds to 6 = c = 2, the two other cases {b,c) = (1,4) or 
(4, 1) were shown in Ref. [10] to be related to the following non-commutative system for R„ € A, 

ff,-.^ Tt (Tt ^-1r pj _/(Rn)'' + l if nodd 



(RnY + 1 if n even 



Let us fix 6 = 1, c = 4. 



Theorem 6.1. The solution R„ of the system (|6.ip is a positive Laurent polynomial of any initial 
data (Rj,Rj_|_i), for any i,n gTj. 

We sketch the proof below. As in the Ai case, one finds two conserved quantities and a linear 
recursion relation for R„ with constant coefficients. Finally, the generating function for R2nR(7^ 
reads: 

Fo(i) = R2„Ro ' = (1 - - t^l - tz-sy^z^y^ 

where 



n>0 



Zl = uiuo 1 = (R^i + l)Ro iRiRo 1 



Z2 = (K-zi)zi-C= (R^iRgR^V(l + Rr^)Ro^(l + Rr^))l^^^I^iR-^^ 

— 1"d2 I /-I I "D— 1\"D— 2 



Z3 = K - Zl = Rr Ro + (1 + Rr )Ro ■ 

and the Laurent positivity follows for (Ro,Ri). The mutated expressions in terms of (Rj,Rj+i) 
are again obtained by manifestly positive non-commutative rearrangements of Fo(i), and the 
general Laurent positivity follows from symmetries of the system. 

6.2. Some rank 2A; + 1 cases. For any fixed k € 7^+, we consider the "rank 2k + 1" recursion 
relation for n € Z, for variables u„ € A: 

(6.2) U2„+2fc+lU2n = U2n+lU2„+2A; + 1 

(6.3) U2„+lU2„+2fc+2 = U2„+2fc+lU2n+2 + 1 

Admissible initial data for this system are Xj = (uj, Uj+i, . . . , Uj_|_2fc) for all i S Z. We have the 
following positivity theorem. 

Theorem 6.2. The general solution u„ to the system (I6.2p - (16.3p is a positive Laurent polynomial 
of any initial data Xj, for all n,i ^"L. 

This is proved in the Appendix [Bl by using the integrability of the system and constructing 
non- commutative weighted path models for tt„. The commutative version of the system (I6.2p - (l6.3p 
corresponds to the affine A2k cluster algebra [1]. Note that the exchange matrix of the fundamental 
seed is singular, so there is no quantum cluster algebra version of the system (I6.2|) - (l6.3p . 



6.3. A non-coprime rank 2 case related to rank 3. The following non-commutative recursion 
relation also satisfies the positive Laurent property: 

(6.4) R„+l(R„) """Rn-l = Rn + 1 

Note that we still have the quasi-commutation relation R„R„_(_i = R„-(_iCR„ from the conserved 
quantity C, and that we may equivalently write R„+iCR„_i = R„(R„ + 1), or 

(6.5) R„_i(R„) ^Rn+i = R„ + C ^ 
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In the commuting case, it would correspond to the recursion relation Rn-i-iRn-i = Rn{Rn + 1) 
which does not correspond to a cluster algebra mutation, as the two terms on the r.h.s. are not 
coprime. However, this relation may be connected to a non-commutative rank 3 system, with the 
positive Laurent property, as follows. Introduce the change of variables: 

R'2n+l = U2n+lU2n, R.2n = U2n-lU2„ 

for some variables u„, n G Z. We have: 

U2n+lU2n-2 = '^2n-1^2n + 1 
U2n-lU2n+2 = U2„+iU2„ + C ^ 

This is nothing but the k = 1 system (|6.2p - (|6.3p with an additional fixed coefficient C~^. Choosing 
the initial data uq = Ri, ui = 1 and U2 = R2, with = uoU2^Uq^U2, we may repeat the 
proof of Appendix |B] so as to include the fixed coefficient C~^, and deduce the Laurent positivity 
property in terms of (Ri,R2). 

7. Towards a non-commutative Lindstrom-Gessel-Viennot Theorem 

All the cases we have been investigating so far are expressible in terms of non-commutative 
weighted path models. More precisely, we have been able to represent the generating function 
Sn>o ^"I^nRo ^ the partition function for paths with non-commutative weights. In the ap- 
plication of the Ar case to the quantum Q-system of Section [5l we also have variables Ra,n) 
n > 1, expressed as "quantum Wronskians" (see Remark l5.6p . This is actually a non-commutative 
generalization of the "Wronskian" expressions for the Q-system solutions in the commuting case 
(|'2.'2p . which take the form: 

In the cases m = mg or mi for instance (fundamental initial data and Stieltjes point, both with 
mi = 0) the Lindstrom-Gessel-Viennot theorem ^25j ^20j allows to interpret directly the l.h.s. as 
the partition function for a family of a non-intersecting lattice paths, namely paths traveled by 
walkers starting at the root of the graph at times 0, 2, 2a — 4, 2a — 2 and returning at times 
2n -|- 2a — 2, 2n -|- 2a — 4, 2n -|- 2, 2n, and such that no two walkers meet at any vertex of T^. 
In the case of general m, we have a more subtle application of the theorem, in which the paths 
must be strongly non-intersecting (see Ref. [7] for details). 

It would be extremely interesting to extend the classical Lindstrom-Gessel-Viennot theorem to 
the case of paths with non-commutative weights. Such a generalization was introduced in ^14j 
in the context of non-commutative Schur functions, but does not seem to apply directly to the 
present framework. 

So far, we have found such an interpretation only in the fully non-commutative Ai Q-system 
case. We wish to interpret the Q-system relation R„_(.iCR„_i — R^ = 1 in terms of pairs of 
non-intersecting paths on the graph r(o), i.e. the integer segment [1,4]. Let P„ = R„Rq^ denote 
the partition function of paths from vertex to itself on F^g) with the non-commutative weights 
yi ) y2) ys respectively for the steps 2 — )• 1, 3 — >■ 2 and 4 — >■ 3, while the up steps 1 — )■ 2, 2 ^ 3 and 
3 — )• 4 receive the trivial weight 1. Multiplying the Q-system relation by C, we have: 

(7.1) P„+iRoCP„_iRoC — P„RoPnRoC = C 

To proceed, we note the following important relations between the weights yi,y2,y3 and Rq, C: 

(7.2) (y3)"y2yiRoC(yi)'^ = y2yiRoC = Ro 1 (yi)"Ro(y3)" = Ro (n e Z) 
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Figure 1. A typical contribution to P„+iRoCP„_iRoC for n = 8 is made of a pair of 
paths A ^ B (blue, left to right) and C ^ D (red, right to left), with two intermediate 
stations at vertex 1, — > C and — > A. When paths intersect (like here), this may also 
be interpreted as a contribution to P„RoPnRoC, upon reversing the direction of travel 
of the portions of paths to the right of their rightmost intersection X, with a weight Rq 
for the step C ^ B. 

We now interpret the quantity P„+iRoCP„_iRoC as the partition function for a path loop 
A^B^C^D^A (see Fig. [T]), where the path A ^ B (partition function Pn+i) is traveled 
from left to right and assigned the usual descending step weights for a step i + 1 ^ z, while the 
portions B ^ C , D ^ A receive the weight RqC, and the path C ^ D (partition function Pn-i) 
is traveled from right to left, and receives the usual descending step weights for a step i + 1 — )• i. 
The contribution for the configuration in the top of FigHJis y3y2yiy3y2yiR'oCyfy3y2yiRoC. 

Theorem 7.1. The quantity 

(7.3) Pn+lR-oCPn-lRoC — P„RoPnR-oC 

is the partition function of loops A^B^C^D^A formed by pairs of non-intersecting 
paths on [1,4] respectively from A = (0, 1) to B = (2n + 2, 1), then from C = (2n, 1) to D = (2, 1), 
with the two intermediate stations B ^ C and D ^ A at the root receiving the weight RqC. 

Proof. Recall first that the portions of loop from B ^ C , and D ^ A receive the weight RqC. 
We define the weight for the same portions, but traveled in the other direction [C ^ B and 
A ^ D) io be Rq. Let us assume the two paths B ^ C, and D ^ A intersect, namely that they 
have at least a common vertex. Let us denote by X the rightmost such common vertex. 

We may now decompose the above loop into three pieces: that from A—^X,X^B^C—>X 
and that from X ^ A. The central portion X ^ B ^ C X, with weight Wxbcx, niay be 
"flipped" into X — > C — >■ B — )• X, with new weight WxcBX, by reversing the direction of travel, 
so that the new configuration now corresponds to pairs of paths ^ — > C and B ^ D oi equal 
length 2n, now separated by the two intermediate stations C ^ B and D ^ A with respective 
weights Rq and RqC, thus contributing to P^RoPriH-oC. 

We now note that WxcBX = Wxbcx- Indeed, depending on the height of X, we have the two 
following possibilities (the first case is illustrated in FigH]): 

• ii) X = {2t + 1, 2) for some t, then the blue path X^5is3^3^4^3^4^ > 

4 ^ 3 -> 2 (weight y3'y2yi, m = n-t) and the red path C— ^■Xisl— )'2->1^^2-> 
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■■■^1^2 (weight y™). The weight of the central loop X ^ B ^ C ^ X is: 

WxBcx = yr(y2yiRoC)yr = R^' = yiRoy2 = yr+'Roy™y2 = Wxcbx 

by use of ([721) • 

• (ii) X = {2t + 2, 3) for some t, then the blue path X-;>Sis3-;>4^3^4-> > A ^ 

3 2 (weight y"y2yi, m = n-t) and the red path C->Xisl^^2^1->2->----> 
1 ^ 2 ^ 3 (weight yf ~^). The weight of the central loop X B ^ C ^ X is: 

WxBcx = yr(y2yiRoC)yr-' = ysR^' = Rr' = y2Ro = y2yrRoyr = Wxcbx 
by use of ^Llh . 

The above flipping procedure is a bijection between pairs of paths A B and C ^ D that 
intersect and pairs of paths A ^ C and B ^ D. They cancel each other in the formula (j7.3p . 
leaving us only with non-intersecting paths, and the theorem follows. □ 

Note that there is a unique configuration of non-intersecting paths, where A^B is 1^2^ 
3^>4^3^>---^>4^3^>2^1andC^L)isl^2^1^2^>---^2^1. The 
corresponding weight is y3~^(y2yiRoC)y"~^RoC = C by (|7.2p . and we recover (j7.ip . 

8. Conclusion: Towards non-commutative cluster algebras 

The rank 2 non-commutative cluster algebras are defined as follows. Let us consider a bipartite 
chain T2 (Cayley tree of degree 2), with alternating black and white vertices indexed by t; E Z 
(black if V is even, white otherwise), and edges with alternating labels 1,2, and a fixed element 
C G A. To each vertex v of T2, we attach a cluster, i.e. a vector of the form x^, = (a, b) in A'^, 
where the cluster variables a, b obey the quasi-commutation relations 

ab = bCa if vis black 
b a = aCb if v is white 

To each (doubly oriented) edge labeled i connecting a black and a white vertex, we associate 
the two following mutations: (i) the forward mutation fif corresponding to the oriented edge 
from the black vertex to the white (ii) the backward mutation = corresponding to the 

oriented edge from the white vertex to the black. Depending on the color of the vertex v, the 
forward and backward mutations of = (a, b) are defined as follows: 



black : fif 


:(a 


b) 


I-)- 


(a' 


,b) 


with 


a'Ca = 


(b)'' + 1 




:(a 


b) 


i-> 


(a. 


b') 


with 


bCb' = 


(a)^ + 1 


t> white : /i^ 


:(a 


b) 


i-> 


(a. 


b') 


with 


b'Cb = 


(a)- + l 




:(a 


b) 


i-> 


(a' 


,b) 


with 


aCa' = 


(b)'' + 1 



These definitions are compatible with an exchange matrix equal to B, = ^ J for black 

vertices and Bo = —B, for white ones, evolving according to the rules of the usual rank 2 cluster 
algebras. The general conjecture of [23] may be rephrased as follows: for any positive integers 
b, c, the cluster variables at any node of T2 may be expressed as positive Laurent polynomials 
of the cluster variables at any other node 0/T2. This is easily proved in the "finite cases" for 
be < 4 (leading to only finitely many distinct cluster variables in the commuting case), where one 
can show explicitly that there exists a positive integer m^^c such that x^+m^ = Cx„C~^, with 
mi^i = 5, mi^2 = 6, and m-i_3 = 8, leaving us with only finitely many cluster variables to inspect. 
In the affine case be = 4, this was proved in |T0] (see the Ai case 6 = c = 2 in the introduction, 
and the case 6 = l,c = 4 Theorem 16. II above), while the case 6c > 4 is still open. 
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To generalize this construction, we may think of using the results of the present paper in 
the following way. Assuming that the Ar Q-system cluster algebra introduced in [21] may be 
generalized to a non-commutative setting, we would expect the relevant cluster variables to still 
be expressible in terms of non-commutatively weighted paths on target graphs, in terms of given 
initial data. In this paper, we have shown how to introduce mutations of non-commutative weights 
via non-commutative continued fractions. Concentrating first on the Stieltjes point, we may trace 
the expressions for Ai„ — R-i^n back to the initial variables Aq,q, and Aqq,_i, a — l,2,...,r + 1 
(see Theorem 13. 17p . This seems to indicate, in this case at least, that the quasi- Wronskians are 
good non-commutative cluster variables. However, the general expressions of Theorem 13.331 for 
Ri^n in terms of quasi- Wronskians or the more compact C,D variables (see ()3.39p and ()3.4ip ) do 
not clearly display which initial data they are related to. We clearly need another definition of 
the non-commutative cluster variables. 

We hope to report on this issue in a later publication. 



Appendix A. Relation of weighted graphs Gm to the graphs Fm of [6] 

In p], we introduced a family of weighted graphs constructed from a Motzkin path m. Here 
we recall the definition of these graphs, as well as the sequence of bijections which brings them to 
form of the weighted graphs Gm which we use in the main body of the paper. It is important to 
understand this connection, because the most basic proof of positivity relies on the fact that the 
graphs Tin are positively weighted with weights which are monomials in the initial data R{m). 

In fact, we describe a sequence of simple bijections 

m o Tm ^ r'^ o Gm- 

The graph is the "compactified" graph of Note that in general, only the graph Fm has 
manifestly positive weights. 

Below, let m be a fixed Motzkin path of length r. 



A.l. The graphs T^. 

(1) Decompose the Motzkin path m into a disjoint union of maximal "descending" segments, 
including those of length one, such that m = U(7nj, mj+i, mj+fc._i). Here, mi+a = 

i 

rrii — a, < a < k in the segment beginning with index i and having length k. 

Example A.l. If m = (2, 1, 2, 2, 3, 2, 1, 1) then m = (2, 1) U (2) U (2) U (3, 2, 1) U (1). 

For notational purposes, we call the maximal descending segment starting with rrij m^*^ 
so that m = Ujm*^*\ 

(2) To each maximal descending segment m'^*) of length, say, k, is associated a graph Fjj^{i). 
It has 2{k + 1) vertices, labeled (i, i + 1, i + k) and («', {i + 1)', {i + k)'). It has the 
following list of oriented edges: (a, a + 1), (a -|- 1, a), (a, a'), (a', a) (for all a), and (6, a) for 
allA; + i>6>a + l >2. Below we give a picture of the graph Fj^{i) for the example of 
the descending subpath m^^) = (3,2, 1,0). 



NONCOMMUTATIVE INTEGRABILITY, PATHS AND QUASI-DETERMINANTS 



39 




(In this picture, an unoriented edge on the graph F is considered to be doubly-oriented, 
for ease of reading.) 

(3) Consecutive maximal descending segments are separated by either a "flat" step of the 
form TUa = rua+i, or an "ascending" step, rua + 1 = ma+i- 

Example A. 2. For the Motzkin path in Example \A.1[ the step separating va.^^'^ = (2, 1) 
and m^^) = (2) is ascending, and the step separating m^^^ = (2) from m^'^) = (2) is flat. 

We "glue" consecutive graphs Fi = and F2 = '^rnO) where m*^'-' and wS-^'^ are 

neighboring descending segments in m with j > i. The gluing procedure is dictated by 
whether the segments are separated by an ascending or flat step. We assume m^*-* has 
length ki and m^-') has length kj. 

(a) If the separation between the segments is flat, we identify the edge and vertices 
ifj) € F2 with {{k,+iy,ki+i) G Fi. 

(b) If the separation between the segments is "ascending", we identify the vertices and 
edges G F2 with {{ki + i)', ki + i) G Fi (i.e., the identification is in the opposite 
sense). 

We define Fm to be the result of gluing the graphs corresponding to consecutive strictly descending 
pieces of m according to the procedure given in (4) above. 

Example A. 3. The purely ascending Motzkin path m = mi = (0, 1, A:) maps to the chain with 
2k + 2 vertices, connected by doubly- oriented edges. The flat Motzkin path mo = (0, 0) of length 
k maps to the graph 

{k + iy • • k + 1 



mo — 




As an example of the gluing procedure, in the case of r = 3 and m = (0, 1,0), we have 
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Fm has a root vertex, which is the vertex 1' in the picture of , or the bottom-most vertex 
in Fm in general. 

The weights attached to the graph Fm are as fohows. Each edge connecting neighboring 
vertices and pointing towards the root vertex 1' carries weight ti/a, with the edges numbered 
a = 1, 2, 2r + 1 from the root vertex on (with edge (i', i) preceding (i + 1, i) in cases where it 
exists). The weights are called the skeleton weights of T^. 

Let Wa,b be the weight corresponding to the edge connecting vertex a to b. Then Wa,a+i = 
Wa,a' = 1, Wa+i,a = ty2a and Wa' ,a = ty2a-i- The Weights Wa,b of the edges (a, 6) with {a - b > I) 
are called redundant weights, and are expressed solely in terms of skeleton weights between the 
vertices a, b as: 

a a— 1 

(A.l) Wa,b = t Yl y2i^2 n 2/2~i-i'(a-^ > !)■ 

i = b+l 4 = 6+1 

These graphs have manifestly positive weights. Let m be a Motzkin path with first component 
mi. Then: 

Theorem A. 4. [6] The ratio Ri^n+mi/Ri,mi is the coefficient of in the partition function of 
paths on Fm from the vertex 1' to itself. 

This theorem can be stated in terms of the resolvent of the transfer matrix associated with 
the graph in the following way. Let be the matrix with entries {Tr^)a,b = ^^a,6- Then the 
generating function for paths from vertex 1' to itself on is the rcsolvGiit 

(A.2) F^{t) = ((I-Tr^)-!)^,^^, = 

n>0 

A.2. Compactified graphs. We now define the graphs F^, which have the property that 

i{^-TrJ~')i,i = {{l-TrJ-\i. 

Note that Zy^y{t) = 1 + tyi{{I - Tr'^)"^)i,i. 

The compactified graphs F^ are obtained from Fm as explained in Section 6 of [8j. This 
process involves an identification of pairs of vertices in the graph Fm, such that the final number 
of vertices is reduced to r + 1. This produces loops and also, in certain cases, extra directed edges. 

(1) For subgraphs arising from descending or fiat segments of m, compactified graphs are 
obtained simply from identification of the vertices i and i' for each i, resulting in a loop 
at vertex i with weight Wj'^j. All other edges retain their weights. 

(2) Consider the vertices and edges arising from each ascending segment of length k in m. 
This is a chain of 2 A; + 2 vertices. 
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(a) Name the vertices from bottom to top, 1, 2k + 2 and identify the vertices 2a — 1, 2a 
(a = 1, k + 1). This leaves k + 1 vertices. 

(b) Rename the vertices from bottom to top again, 1, + 1. 

(c) Identification of vertices results in a loop at each vertex a with weight inherited from 
the weight associated with the edge connecting 2a to 2a — 1, W2a,2a-i- The spine 
retains weights W2a+i,2a along each edge (a + 1, a). 

(d) In addition we adjoin ascending directed edges, connecting vertex a to 6 for all 1 < 
a < b < k + 1, with weight (—11 



\b-a+l 



Remark A. 5. In fact, the construction of the compactified graph directly from the Motzkin path 
is easy to describe. The Motzkin path is decomposed into maximal ascending, descending and 
flat segments, the corresponding compactified graphs and their weights are constructed as above, 
and the graphs are glued: The top two vertices corresponding to a segment are identified with the 
bottom two vertices corresponding to the next segment. Incident edges between these vertices are 
identified. (See pictures below.) 

Lemma A. 6. [8] We have the equality of partition functions 

{i^-TrJ~'),^, = {{I-TrJ-'),, 

In light of remark [A. 51 we can focus on graph segments corresponding to segments of m which 
do not change direction. 

A. 3. The graph Gm- We now turn to the equivalence of resolvents of the graph Gm, defined in 
Section 2, to those on the graph T^. The difference is the presence of long directed edges in 
(connecting vertices along the spine of the graph with indices which differ by 2 or more) . 

On the level of the transfer matrix Tp' , these edges are represented by non-zero elements above 
or below the second diagonals. These entries can be eliminated using row (column) elimination 
for descending (ascending) edges, respectively. Conveniently, 

Lemma A. 7. The matrices E and F such that E{1 — T-pi^)F has no non-trivial entries away 
from the second diagonals are such that E is upper unitriangular and F is lower unitriangular. 



The (1, l)-resolvent 
is obviously unchanged under such transformations. 



Proof. It is sufficient to show this for each block of T' = Tp' corresponding to maximal descending 
and maximal ascending pieces of m. The case of flat pieces is trivial as it has no entries off the 
second diagonals. 

(1) Descending pieces. The matrix I — T' related to the vertices in the graph corre- 
sponding to a descending path (k — l,k — 2, 0) is the {k + 1) x (k + l)-matrix 



/ l-xi 



I-T' 



-X2 
X2X4 
X3 

2:3X5 



1 



-1 

- 

— X4 

X4,Xg 

xs 





-1 

- X5 

-xe 







\ 





-1 

1 — Xj 

-X2k-2 

X2k-2X2k 
X2k-1 



X2k-1 
-X2k 



-1 



X2k+1 / 



42 



PHILIPPE DI FRANCESCO AND RINAT KEDEM 



Here, X2a-i = Wa'^a and X2a = w^a+i,a- Clearly, there is an lower triangular matrix which 
eliminates all of the entries below the second diagonal when multiplied on the right: 

X2i 



1,(1 <i <k + l) 



The resulting matrix M 
i<k-l. 



Fj+i i = —{1 < k — 1), Fi j = otherwise. 

X2i+1 

(I — T')F has non-vanishing entries Mj^j+i = — 



1, and for 



(A.3) 



1 - {X2i-1 



X2i 



M, 



X2i 



X2i+1 X2i+l 

In addition, the last two columns of the matrix M are those of I — T'. 

Note that we have lost the positivity property of the weights in this case, although the 
generating function of paths is positive. Note also that the weights attached to the edges 
connecting the last two vertices have not changed under this transformation. 
(2) Ascending pieces. The matrix T', attached to the vertices in the graph correspond- 

., /c — 1) is the (k + 1) X {k + l)-matrix 



ing to an ascending path (0 
/ l-xi 

-X2 



I-T' 



-1 


1 -1 


1 


±1 


1 - X3 


-1 1 


-1 




— X4 


1-X5 -1 


1 







— Xq 1 — Xj 


-1 





-X2k-2 i - X2k~l 
\ • • • -X2k 1 - X2k+l J 

There is obviously an upper triangular matrix corresponding to row operations on 
this matrix, or multiplication on the left, such that M = EiJ. — T') has no entries away 
from the second off-diagonals. The resolvent is equal to the resolvent of the matrix with 
non- vanishing entries Mj^j+i = — 1, and if i < k — 1, 

(A.4) Mii = l- {xi + Xj+i), Mi+i^i = -X2iX2i+i. 

The last two rows of M are those of I — T'. Note that in this case, we retain the positivity 
property of the weights. 

□ 

To summarize, we can draw the graphs corresponding to the sequence of transformations from 
m (an ascending, descending or flat path) to in graphical notation: 

(1) For ascending paths. 



1 



-1 



2k + 2 




2k + 1 



rn 




Q 

O 



fc + i 



Gr 
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(2) For descending paths, 




J- m -p/ (~i 

In the final graphs Gm for ascending and descending Motzkm paths, the symbols x and = along 
certain edges indicates that the weights have been renormalized from the initial values X2i-\ and 
X2i to the values X2i-\ — X2i/x2i+i and X2i/x2i+i for = (see Eq. (lA.3|l ) and to the values X2i-i + X2i 
and X2iX2i+i for x (see Eq. (IA.4|) ). For Xi = tin with the appropriate labeling of skeleton weights, 
these are precisely the weights y of Definition 12.11 of Section [2j 

Gluing two path pieces, m followed by m' into a longer Motzkin path corresponds to identifying 
the two top vertices (and incident edges) of Gm, the lower graph, with the bottom two of Gm', 
the top graph. In the cases where the weights have been renormalized in the lower part of the 
upper path, the renormalized weights appear in the glued graph. 



Appendix B. Proof of Theorem 16.21 

Due to the invariances of the system ()6.2p - (j6.3p . we only need to express the solutions u„ in 
terms of the initial data xq = (uo,ui,U2, ...,U2fc). Indeed, introducing the anti-automorphism ip 
such that (/?(!) = 1, <^(a6) = ip{b)ip{a) for all a,b € A, and (/^(uj) = Uj+i, for i = 0,1,. ..,2k, we 
easily find that (^(u„) = u„+i for all n G Z. So if we have an expression u„ = /„(uo,ui, ...,U2k) 
for the general solution of (j6.2p - (j6.3p in terms of the initial data xq, we may apply if iteratively 
on it to get u„ = 99*(/„_i(uo, ui, U2fc)) = g„,j(ui, Uj+i, . . . , Uj+2fc), namely an expression for u„ 
in terms of any other initial data Xj, for all n,i S Z. We may also restrict ourselves to u„ with 
n £ Z-|.. Indeed, let ip be the anti-automorphism such that ■0(uj) = U2k~i-i, i = 0, 1, 2k, then 
we have ip{un) = U2k-i-n for all n S Z. Hence if we have an expression u„ = /„(uo,ui, ...,U2k) 
for n > 2k + 1, then applying iptlt to it gives U2k-n = ^n(uo, ui, U2k), namely an expression for 
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u„, n < 0, in terms of xg. Moreover all of the above expressions are positive Laurent polynomials 
iff /„ is a Laurent polynomial for n > 0. 

By the above remarks, we only need prove the statement of the theorem for n > and i = 0. 
We first need the following lemma, expressing the integrability of the system ()6.2p - ()6.3p . 

Lemma B.l. There exist two linear recursion relations of the form: 

U2n+2fc+l ~ U2„+lK + U2„_2fc+1 = 
U2n+2A: — Ku2„ + U2„_2fc = 

for all n ^"L, and with K expressed in terms of xq as: 

k 

(B.l) K = UQU-^i^ + U2fcUo ^ + ^ U2/_i (u^/ + U^/+2) 

Proof. Introduce 

Kn = ^2n+li^'in+2k+l + U2n-2fc+l) = (U2n+2fc + U2n-2fe)u2„^ 

Subtracting (fST^ for n n-k from that for n, we get (u2„+2fe+i+U2n-2fc+i)u2n = U2„+i(u2„_2A:+ 
U2n+2fc), hence K„ = L„. Moreover, 

U2n+l-f^nU2n+2 = U2n+lU2n+2fc+2 — 1 + U2„-2A;+lU2n+2 

= U2n+l(u2„+2A:+2 + U2n-2A:+2) = '^2n+lLn+l'^2n+2 

hence K„ = L^+i = K„+i = K is a conserved quantity of the system ()6.2p - (|6.3p . Eq. (|B.ip is 
easily derived by induction on fc, and the lemma follows. □ 

Note that (/j(L„) = K„, hence ^{K.) = K. Analogously, ^(L„) = 'Kk-i-n-, hence ^'(K.) = K as 
well. 

Introducing the generating functions 

oo oo 

"^2i+2kn l+2fcn 
ra=0 n=0 

for i = 0, 1, k — \ and the weights: 

X2i = U2j+2fcU2/ 

Z2j = K — X2j = \l2i-2k'^2i 

W2i = Z2iX2i - 1 

X2i+1 = 'il2i+l'^'2i+l+2k 

Z2i+1 = K — X2j+l = U2j^,_;^U2i+i_2/!c 

W2i+1 = X2i4.lZ2i+l — 1 

Lemma IB . 1 1 implies the following continued fraction expressions: 

Fi{t) = {l-{^2^+Z2^)t + t^)-\l-Z2it)u2i 

= (l - X2it - (1 - Z2jt)"V2it^) ^\l2i 

Gi{t) = U2i+l(l - Z2i+lt)(l - (X2i+l + Z2i+l)t + t^)"-^ 

= U2i+l (l - tX2i+l - t^W2i+l(l - Z2i+lt)"-^) ^ 

whose series expansions in t have coefficients that are manifestly positive Laurent polynomials 
of (xj,Zj,Wj). Finally, we see that the x's are defined recursively by xq = U2fcU^^ and X2j = 
X2j-i + U2^^.^_-|^U2j.^, X2j+i = X2j + ^2j^^i^2j^ ' hencc they form subsums of K, and so do the z's. 



NONCOMMUTATIVE INTEGRABILITY, PATHS AND QUASI-DETERMINANTS 



45 



Also, all x's contain the term U2A;Uq , while all z's contain uqu^^ , hence both XjZj and ZjXj contain 
1: we conclude that the x, z,w's are all positive Laurent polynomials of xq. This completes the 
proof of Theorem 16.21 

The existence of the linear recursion relations of Lemma IB. II implies by Proposition 13.131 that 
all 3 X 3 quasi- Wronskians '^2kn+i — |(u2fcn+i-2a-26+8)i<a,fe<3|i i rnust vanish. Each function 
Fj(t), i = 0, 1, A: — 1, is actually a particular instance of the general non-commutative case 
with r = 1, discussed in Section [3l It is the generating function of paths on the graph Gi, with 
respective weights wi_2 = 1 and: 



As we also have a condition W2j = Z2jX2i — 1, we deduce that the conserved quantity C2,n = 
C = ysyi = 1. This provides us with examples in the general Ai setting that are different from 
the non-commutative Ai Q-system of Section II. 3[ In the latter indeed, C is kept arbitrary, but 
the value of A2,n = H-n+i — RnR-^liRn is fixed by the Q-system relation A2,n = R-^^'^R-n-iRn- 
Here we have C = 1, hence the recursion relation /S.2,n+i = H-n-iR'n^^2,ra (with R„ = U2i+2fcn)) 
but there is no such additional relation. Note that C = 1 also implies that there is no quantum 
version in which C would be a non-trivial central element of A. This is in agreement with the 
above-mentioned fact that the exchange matrix of the associated cluster algebra is singular. 
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